Phase coherent transmission through 
interacting mesoscopic systems 

Gregor Hackenbroich 
Universitdt GH Essen, 45117 Essen, Germany 



Abstract 

This is a review of the phase coherent transmission through interacting 
mesoscopic conductors. As a paradigm we study the transmission amplitude 
and the dephasing rate for electron transport through a quantum dot in the 
Coulomb blockade regime. We summarize experimental and theoretical work 
devoted to the phase of the transmission amplitude. It is shown that the evo- 
lution of the transmission phase may be dominated by non-universal features 
in the short-time dynamics of the quantum dot. The controlled dephasing 
in Coulomb coupled conductors is investigated. Examples comprise a single 
or multiple quantum dots in close vicinity to a quantum point contact. The 
current through the quantum point contact "measures" the state of the dots 
and causes dephasing. The dephasing rate is derived using widely different 
theoretical approaches. The Coulomb coupling between mesoscopic conduc- 
tors may prove useful for future work on electron coherence and quantum 
computing. 
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I. INTRODUCTION 



Electron transport on mesoscopic length scales has been intensively studied during the 
last fifteen years |1]-|4]] . A variety of new phenomena has been discovered such as the weak 
localization, universal conductance fluctuations, Aharonov-Bohm oscillations in mesoscopic 
rings, or persistent currents. The origin of these phenomena is the quantum mechanical 
phase coherence of the electronic wave functions. The degree of coherence can be measured 
by the phase coherence length which is the typical length on which electrons travel 
without loosing their phase coherence. In mesoscopic systems, exceeds the system size, 
and transport exhibits quantum interference effects. 

How is mesoscopic transport modified by electron-electron interactions? Interactions 
generally play a minor role in good conductors where electron wave functions are delocalized 
over the whole system and screening is effective. Such systems, both in the presence and 
absence of disorder, have been successfully described by noninteracting electron models 
||. Interactions only lead to small corrections, e.g. they suppress the phase coherence 
of interfering electrons by coupling them to the bath of other electrons. In contrast to 
the transport in weakly interacting conductors, interactions strongly affect transport in 
high magnetic fields and in spatially confined geometries. Important examples are the 
fractional quantum Hall effect and the tunneling through small spatially confined electron 
islands, known as quantum dots. The electron interactions in a quantum dot give rise to 
a characteristic energy scale called the charging energy. Similar to the ionization energy 
of an atom, the charging energy is the energy necessary to remove or add a single electron 
to the quantum dot. The charging energy in state-of-the-art semiconductor quantum dots 
typically exceeds the single-particle level spacing by a factor 10 or more. Charging effects 
can drastically alter and nearly suppress the transport of charge through these systems at 
temperatures roughly below 1 K. 

In this work we review the phase coherent transmission through strongly interacting 
systems. As a paradigm, we study the tunneling through quantum dots. The subject 
started in 1995 with an ingenious experiment by Yacoby et al. 0. The authors utilized a 
novel interference device with a quantum dot embedded in one arm of an Aharonov-Bohm 
(AB) ring (see Fig. |l|). The measurement of flux-periodic current oscillations through the 
device proved for the first time that part of the tunneling current through a quantum dot 
is coherent. In a step beyond the demonstration of coherence, an improved version @ of 
the experiment allowed to measure the phase of the transmission amplitude through the 
quantum dot. At sufficiently low temperatures, this phase is carried by a single resonant 
many-body state. The transmission phase displayed a number of unexpected properties. 
Most notably, virtually the same transmission phase was found for a whole sequence of 
resonant quantum dot states; between neighboring states the phase displayed a sharp phase 
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FIG. 1. Electron micrograph of the device used in the experiment ||. The brighter regions 
indicate metallic gates; B is a metallic air bridge. Electrons flow between the regions labeled S and 
D through the left or the right arm of an Aharonov-Bohm ring. The quantum dot is inserted in 
the left arm. Taken from Ref. ||. 

slip. These features have been addressed but not yet fully explained in a large number of 
theoretical papers. 

The field expanded substantially in 1998 when it was realized that the coherence of quan- 
tum dot states can be controlled by external means. Controlled decoherence was achieved 
in a device with a quantum dot that was capacitevly coupled to a quantum point con- 
tact in close vicinity. The quantum point contact acted as a measuring apparatus for the 
number of electrons on the quantum dot. Since number and phase are conjugate variables, 
the measurement caused the dephasing (decoherence) of electron states in the quantum dot. 
The loss of coherence was detected as a suppression of AB oscillations across the quantum 
dot. The experiments for the first time allowed to study fundamental principles of 
quantum mechanics such as the number-phase or the particle-wave duality with solid state 
devices. New experiments in this direction, e.g. testing the suppression of tunneling due to 



frequent measurements, have been proposed [ 1C , 11 



The present work reviews the main theoretical and experimental results on phase co- 
herent transport through quantum dots. We emphasize the theoretical developments. The 
most important formal tool is scattering theory. We use both the well-known scattering 
approach of Landauer to discuss aspects of noninteracting electrons and many-body scat- 
tering theory to address interaction effects. The scattering methods are supplemented by 
other techniques, most prominently by master equation approaches. Most of the theoretical 
work was motivated by and must be discussed in connection with the beautiful experiments 
I^H on quantum coherence done by researchers at the Weizmann Institute. An account of 
the relevant experimental results, therefore, precedes the theoretical discussions. It is as- 
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sumed that the reader is familiar with ballistic transport in mesoscopic systems. We assume 
basic knowledge of the Landauer formula, and of transport through Aharonov-Bohm rings, 
quantum point contacts and quantum dots. An overview about most of these concepts can 
be found in Refs. p]-§J. For a recent review on transport through quantum dots in the 



Coulomb blockade regime see Ref. Reference |L3[ gives a popular account of quantum 
dots and the Coulomb blockade. 

In Sec. [TI| we summarize the first experiments on the phase coherent transmission 
through quantum dots in the Coulomb blockade regime. Both experiments aimed at mea- 
suring the phase of the transmission amplitude through a quantum dot. It turned out, that 
in the first experiment this measurement was hampered by peculiar symmetries resulting 
from the two-terminal set-up used in that experiment. Section [ill] is devoted to a discussion 
of these symmetries. We derive the two-terminal conductance through an AB ring with a 
quantum dot. The theory explains the sharp jump of the AB oscillation pattern found in the 
vicinity of transmission resonances. Section |IV| deals with the observed phase similarity of 
transmission peaks and with the sharp phase jump between the peaks. We show that these 
observations may be traced back to features of the short-time dynamics in the quantum dot 
that completely dominate the transport in the tunneling regime. In Sec. |V| we discuss the 
controlled dephasing of quantum dot levels due to the coupling to a mesoscopic conductors 
in the environment. Section |V| deals with new theoretical proposals to study the dephasing 
of charge oscillations in coupled quantum dots. The proposed experiments could provide ev- 
idence for a fundamental quantum effect known as the quantum Zeno effect: the suppression 
of quantum transitions due to frequent observations with a measuring device. We conclude 



in Sec. VII. 
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II. MEASURING THE ELECTRON TRANSMISSION PHASE 



The phase coherence of electronic wave functions lies at the heart of mesoscopic interfer- 
ence phenomena like the universal conductance fluctuations or weak-localization. In most 
cases, quantum interference has been studied with mesoscopic conductors that were con- 
nected to few or many conducting channels. Interaction effects are usually small in these 
systems due to good screening of the electron charge. Charging effects become important 
when the conductor is small and only weakly coupled to the external leads. A paradigm is 
a quantum dot in the Coulomb blockade regime. 

The idea and the practical implementation of devices for measuring the electron trans- 
mission phase through a quantum dot was developed in a series of ingenious experiments ||[7j 
at the Weizmann Institute. The experiments utilized an Aharonov-Bohm (AB) ring with a 
quantum dot in one of its arms. The quantum dot was in the Coulomb blockade regime. The 
measurements of the electron phase revealed a number of unexpected results. The technique 
for measuring the electron phase and the main experimental results are described in this 
section.^] 

A. Experiment of Yacoby et al. 

The first experiment addressing the electron phase through a quantum dot was realized 
by Yacoby et al. ||. The experiment utilized a novel device (see Fig. 0) to measure the phase 
evolution through the dot against a fixed reference phase. The quantum dot was inserted in 
one arm of an Aharonov-Bohm ring. The basic idea was to extract the transmission phase 
through the quantum dot from the phase of the Aharonov-Bohm current across the ring. 
The device was defined by metallic gates on the top of a GaAs-AlGaAs heterostructure. 
The quantum dot was placed in the left arm of the ring, its conductance could be adjusted 
by two quantum point contacts. A third gate (the plunger gate P) controlled the area and 
the electrostatic potential at the dot. The dot was about 0.4/xm x 0.5/xm in size. A special 
lithographic process, invoking a metallic air bridge (B), enabled to contact the center metal 
gate that depleted the electrons underneath the ring's center. Each of the arms of the 
Aharonov-Bohm ring supported a few conducting modes. The ring was connected to two 
external contacts, source (S) and drain (D), between which a small voltage was applied. 



1 We do not discuss the experiment [14 that measured the phase and magnitude of the reflection 



coefficient of a quantum dot. This experiment 14] was performed in the integer quantum Hall 
regime and used a second quantum dot as an interferometer. 
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FIG. 2. Electron micrograph of the device used in the Yacoby-experiment. The brighter regions 
indicate the metal gates. Electrons flow between source and drain through the left or the right arm 
of the Aharonov-Bohm ring. The quantum dot is inserted in the left arm. The central metallic 
island is biased via an air bridge (B). Taken from Ref. ||. 

The elastic mean free path in the two-dimensional electron gas was about 10 fim while 
the diameter of the Aharonov-Bohm ring was ~ 1 — 1.5 /im. The phase coherence length 
was larger than the ring's circumference. The quantum dot had a resistance of 1MQ and a 
very small capacitance C ~ 160 aF corresponding to the charging energy e 2 /2C ~ 0.5 meV. 
The dot contained around 200 electrons. Its average single-particle level spacing was A « 40 
fieV. The experiment was performed at a temperature around 100 mK corresponding to 
the thermal energy kT pa 9 [ieV . The intrinsic width T of the Coulomb peaks was estimated 
from the conductance peak height to T ~ 0.2 fieV. These scales imply that the quantum 
dot was in the Coulomb blockade regime, and that the transmission at each Coulomb peak 
resulted from resonant tunneling through a single or few levels of the quantum dot. 

The experiment |J demonstrated for the first time that part of the tunneling current 
through a quantum dot is coherent. The experimental evidence is depicted in Fig. [| Shown 
is the ring current vs. plunger voltage Vp for a fixed small source-drain voltage. The ring 
current was obtained by subtracting from the total current across the ring a large Vp- 
independent background due to the right arm. The Coulomb blockade in the dot creates 
sharp conductance peaks in the ring current for fixed magnetic field. Fixing the voltage Vp 
on the side of a current peak and sweeping the magnetic field, one observes periodic current 
oscillations. The period of the oscillations corresponds to one flux quantum threading the 
area of the ring, in agreement with the expected Aharonov-Bohm period. The oscillation 
contrast, denned as peak-to-peak current over the average current, was in the range 0.2 to 
0.4. 

In a next step, the current oscillations were investigated at different values of the plunger 
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FIG. 3. One of the ring's current peaks as a function of the plunger gate voltage. The large 

current of the right arm has been subtracted. The top left part shows the Aharonov-Bohm oscil- 
lations of the current vs. magnetic field at fixed Vp = V m . Inset: Oscillation contrast defined as 
peak-to-peak vs. average current through the dot. Taken from Ref. [||. 



voltage Vp. A change in Vp was expected to modify both the magnitude and phase of the 
transmission amplitude through the dot. The experiment was motivated by the idea that the 
change in the transmission phase would be directly reflected in a shift of the Aharonov-Bohm 
oscillations which, in turn, could be seen experimentally. This one-to-one correspondence 
is suggested by the following argument: Suppose the ring and the leads support only one 
conducting channel. According to the Landauer formula, the zero-temperature current 
between the leads is proportional to the ring transmission coefficient \t(Ep)\ 2 at the Fermi 
energy E F . For fully coherent transport, t = t R exp(27rz$/$ ) + t L , where t R and are the 
transmission amplitudes through the right and left arm, respectively, $ is the flux through 
the ring and <3>o the flux quantum. This yields the interference term 



2Re{t L ^exp(-27T2$/$ )} = 2|t fl ||t L | cos(£ L - £r - 27r$/$ ), 



(2.1) 



where £r = arg(t^), £l = arg(t^). Any shift in the phase £l — £r should thus be directly 
reflected in a shift of the Aharonov-Bohm oscillations. 

The above argument turned out to be incorrect. It neglects multiple reflections through 
the ring. The argument fails in particular for a two-terminal geometry, i.e. a ring connected 
to two external leads. It was realized pTOj] shortly after the Yacoby-experiment that 
Onsager symmetries valid for a two-teminal device restrict the phase of the Aharonov-Bohm 
oscillations to either or 7T, spoiling the correspondence between the Aharonov-Bohm phase 
and the transmission phase through the quantum dot. A discussion of this issue is presented 
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FIG. 4. (a) A series of three Coulomb peaks and (b) the current oscillations measured at the 
marked points A, B, and C. All oscillations are seen to be in phase. The large current of the right 
arm has been subtracted. Taken from Ref. ||. 

in Sec. |T|. Despite of its failure for a two-terminal ring, the simple phase-argument catches 
the essential idea for measuring the electron phase as it was later realized in the multiple- 
terminal device of Schuster et al. [|7| . 

Figure |] shows the ring current and the Aharonov-Bohm oscillations measured at three 
successive peaks. The oscillations at similar points (A, B, and C in the figure) have the same 
Aharonov-Bohm phase. This repetition of the phase was found within the whole sequence 
of Coulomb peaks (comprising 12 peaks) investigated in the experiment. The evolution 
of the Aharonov-Bohm phase along a single Coulomb peak is displayed in Fig. |5|. Four 
different interference patterns taken at the points 1, 2, 3, and 4 specified in Fig. ||](a) are 
shown in Fig. |5](b). The Aharonov-Bohm phase of the patterns shifts by 7r as one sweeps 
through the peak. The shift happens rather abruptly between the points 2 and 3. This 
can be seen in Fig. |5](c) which summarizes the phase measurements along a Coulomb peak. 
The sharp jump in the measured Aharonov-Bohm phase is contrasted with the expected 
phase evolution of the transmission amplitude assuming resonant tunneling through a single 
level of the quantum dot. The latter phase increases smoothly on the scale of the peak 
width (which is of order kT). Theoretical arguments [|T^-|TP[ proved that there is no scale 
associated with the rather sharp jump seen in the experiment. 



B. Experiment of Schuster et al. 

The Aharonov-Bohm phase in a two-terminal device is restricted to or it. No such 
restriction exists in a multi-terminal probe. This motivated Schuster et al. [|7| to perform an 
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FIG. 5. Evolution of the Aharonov-Bohm phase along a current peak, (a) Current vs. gate 
voltage at a current peak, (b) A series of interference patterns taken at the points specified in (a). 
The phase jumps between patterns 2 and 3. (c) Phase measured at two peaks (circles and triangles). 
The broken line is guide to the eye. The expected behavior of the quantum dot transmission phase 
in a ID resonant tunneling model is shown by the solid line. Taken from Ref. ||. 

interference experiment similar to the Yaco by-experiment || but with more than two leads 
connected to the interferometer. The electron micrograph of the device and a schematic 
description of the experiment are shown in Figs. |], [?[ The central element of the device is 
an Aharonov-Bohm ring with a quantum dot embedded in its right arm. Several contacts are 
connected to ring, namely the emitter (E), the collector (C) and a base region (B). The base 
contacts were held at zero chemical potential. Incorporated in the structure are additional 
barriers. They reflect diverging electrons into the two-slit device and towards the collector. 
The reflectors were necessary to enhance the collector signal that could otherwise not be 
measured. All contacts were defined by negatively biased gates on top of the heterostructure. 

The quantum dot contained roughly 200 electrons with a mean single-particle level 
spacing around 55 [ieV . The temperature of the two-dimensional electron gas was 100 mK. 
The intrinsic (zero-temperature) width T of the Coulomb peaks was estimated to be of the 
order or even larger than kT. This was achieved by slightly opening the point contacts 
between the quantum dot and the ring. The large value of T corresponds to a quantum 
dot resistance larger than the resistance quantum Rk = h/e 2 but much smaller than the 
dot resistance in the Yacoby-experiment. Working with a modest quantum dot resistance 
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FIG. 6. Scanning electron micrograph of the double-slit device used in the Schuster-experiment 
0. The grey areas are metallic gates on the top of the heterostructure. The quantum dot is inserted 
in the right slit. Taken from Ref. 0. 




FIG. 7. Schematic of the device structure of the Schuster-experiment. An Aharonov-Bohm 
ring is coupled to an emitter (E), a collector (C) and a base region (B). Reflector gates reflect 
diverging electrons towards the collector. The quantum dot is defined by the central electrode and 
the three electrodes on the right hand side of the structure. Taken from Ref. @. 
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enabled Schuster et al. to measure the Aharonov-Bohm oscillations not only at the Coulomb 
peaks but also in between. 

Schuster el al. investigated the voltage drop Vcb between collector and base for a fixed 
excitation voltage Veb applied between the emitter and the base. The Aharonov-Bohm 
interference in the transmission coefficient Tec leads to an oscillatory contribution to Vec 
from which the Aharonov-Bohm phase is extracted. Fig. |5](a) shows Vec measured as 
a function of the gate voltage Vp for a fixed magnetic field. One observes pronounced 
resonance peaks as expected for a quantum dot in the Coulomb blockade regime. When 
the magnetic field is changed, the collector signal shows AB oscillations with the expected 
period AB = §o/A where A is the area of the AB ring. The observed oscillation patterns 
measured at the four points 1, 2, 3, and 4 close to a resonance are shown in Fig. ||(b). The 
oscillation pattern shifts smoothly as one moves through the resonance. Fig. |8|(c) displays 
the phase and the squared magnitude of the AB signal at a resonance peak. The data points 
are represented by full circles. The phase shows the expected monotonic rise by 7r over the 
width of the resonance. 

Schuster et al. compared their data with a theoretical model for resonant transmission. 
They described the coherent part £qd of the transmission amplitude through the dot by the 
Breit-Wigner ansatz 

tQD = iCN E F -E„ + iT N /2' (2 ' 2) 
where Cat is a complex amplitude, Ep the energy of the electrons transmitted through the 
device, and E N and T N the energy and the width of the resonance in the quantum dot. 
Both the squared magnitude |£qd| 2 and the phase argtQD are compared with experimental 
data in Fig. |8](c). Good agreement is found when « 4kT is used as a fit parameter. 

The collector voltage Vcb, the magnitude, and the phase of the AB oscillations measured 
over a sequence of five peaks are shown in Fig. [| The striking observation is that the phase 
is very similar at all peaks. The phase increases roughly by tt at each peak. Note that the 
peaks widen and start to overlap as the plunger voltage increases. At the same time, the 
overall variation of the phase is reduced. The likely origin of these effects is the electrostatic 
influence of the plunger on the point contacts at the quantum dot. They open slightly with 
increasing plunger voltage. A striking feature of the data is the sharp phase lapse by it 
between the resonances. The phase lapse occurs when the magnitude of the AB oscillations 
vanishes. Unlike the phase jump at resonance found in the Yacoby-experiment, there is a 
scale associated with the phase lapse between the peaks. This becomes evident for increasing 
Vp when the lapse smoothens and is resolved experimentally. While Schuster et al. were 
able to model the phase evolution at resonance, they could not explain the phase behavior 
between resonances. The origin of the phase lapse, the associated energy scale, and the 
similarity of the phase at subsequent peaks is addressed in Sec. [TV. 
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FIG. 8. Conductance and phase evolution along a Coulomb peak, (a) Resonance peaks as 
a function of the plunger gate voltage, (b) A series of interference patterns taken at the points 
specified in (a), (c) Squared magnitude and phase of the Aharonov-Bohm oscillations (dots). The 
solid and dashed line are fits for the phase and the squared magnitude obtained with a Breit-Wigner 
model. Taken from Ref. jjj. 
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FIG. 9. (a) A series of Coulomb peaks; (b) Magnitude of the Aharonov-Bohm oscillations; (c) 
Phase of the Aharonov-Bohm oscillations as a function of plunger gate voltage. The solid lines are 
guides to the eye. Taken from Ref. M. 



III. PHASE EVOLUTION IN A TWO TERMINAL DEVICE 



Electrons that tunnel through a quantum dot interact with many other particles including 
electrons in the dot and in the surrounding gates, phonons in the substrate, photons from 
the environment, etc. The interactions may spoil or even destroy the phase coherence of the 
transmitted electrons. The Yacoby-experiment demonstrated for the first time that part 
of the tunnel current is coherent. As discussed in Sec. ||, the experiment employed a quantum 
dot in an AB ring. Phase coherence was demonstrated by measuring AB oscillations of the 
ring current. The oscillations displayed a rigid phase, the phase taking on only two values 
and 7r with abrupt jumps between these values. In this section, we review theoretical 
explanations |To}-|H|] of this observation. We show that the phase rigidity is imposed by the 
two-terminal character of the Yacoby-experiment (Sec. A| ) . Information about about the 
transmission phase through the quantum dot may still be obtained from the amplitude of 
the AB oscillations (Sec. pTl~B| ). 
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A. Phase rigidity 



The argument for the phase rigidity was presented by Levy Yeyati and Biittiker 
Their argument is based on reciprocity relations for the transmission amplitude derived by 



Biittiker 20,21 1. These relations state that the two-terminal conductance is an even function 



of the flux through the ringj^j 



G($) = £(-$). 



(3.1) 



This relation holds for an arbitrary number of transverse channels in the leads. We briefly 
sketch a derivation of Eq. (|3.1|) , valid for one channel in each lead: Time-reversal invariance 
implies ti 2 (Q) = t 2 i(— $) for the transmission amplitude between the two leads labeled 1 
and 2. The unitarity of the scattering matrix yields |ii2($)| 2 = |i2i( < ^ ) )| 2 - Combining these 
relations, one finds |£i2($)| 2 = |ti2 ( — $)| 2 - Substitution in the Landauer formula yields the 
claimed symmetry G($) = G(— $). 

The conductance through the AB ring is a periodic function of <£>. Thus, G can be 
expanded in a Fourier series 



= G + G 1 cos(2tt$/$ + (Ji) + 



(3.2) 



where the dots indicate the higher harmonics. Combining the expansion ( |3.2|) with the 
symmetry (|3.1| ), one finds that the phase S\ can only be or tt (or equivalently an even or 
odd multiple of tt). In the Yacoby-experiment, all expansion coefficients are a function of 
the plunger voltage Vp at the quantum dot. If the phase 5\ changes upon variation of Vp, 
this change must be a sharp jump by tt. 



The above explanation of the phase rigidity relies on the symmetry relation (|3.1|). This 



relation remains valid even in the presence of inelastic scattering [15] . Inelastic processes can 



be included adopting a model suggested by Biittiker |2TJ: In addition to the two physical 
leads 1 and 2, the model includes a fictitious third lead that connects the AB ring with a 
phase-randomizing reservoir. The condition that no current flows through the third lead 
determines the chemical potential of the attached reservoir. The two-terminal conductance 
is found I^TJ to be 



G 



h 



T 12 + 



23^31 



T 32 +T ; 



31 



(3.3) 



where T^- is the transmission coefficient from lead j to lead i. The term T 12 describes the 
direct elastic transmission between leads 1 and 2 while the second term on the right hand side 



2 This symmetry was, implicitly, formulated in earlier work of Gefen, Imry, and Azbel ]22j 
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FIG. 10. Amplitude of the AB harmonics vs. the phase 6. The quantum dot displays transmis- 
sion peaks around 8 = 70 and 8 = 170. Taken from Ref. [18|. 



of Eq. ( |3.3| ) accounts for the inelastic transmission. Using the conservation of transmission 
probability S n T mn = 1 (for a single channel in each lead), one can write Eq. ( [3.3|) in the 
form 

T13T31 



G 



h 



(3.4) 



I-T33. 

The reciprocity relations T mn (<$>) = T nm (— $) imply that G is an even function of flux. 
Consequently, the phase rigidity persists even in the presence of inelastic scattering. The only 
effect of inelastic processes is to reduce the amplitude of the AB oscillations by decreasing 
the direct elastic transmission Tvi- 

The nature of the AB oscillations close to a phase jump was studied both theoretically 
and experimentally by Yacoby, Schuster and Heiblum flSfl . These authors showed that the 
phase jump occurs at the point where the amplitude of the h/e periodic AB oscillations 
vanishes. At that point, the ring conductance is dominated by h/2e periodic oscillations. 



The theoretical analysis in Ref. [18| is based on a formula derived by Gefen, Imry and Azbel 



22] for the two-terminal conductance across a single-channel AB ring with single scatterers 



in each arm, 



G 



2e 2 
~h 



2(Ae^ + Be'**) 



(3.5) 



De 2i* + Ee- 7 ** + C 

Here <fi = 7r$/$ is the dimensionless flux through the ring, and the constants A,. . . ,E de- 
fined in Ref. |22j are functions of the scattering amplitudes of the scatterers in the ring. The 
quantum dot in one of the arms is modeled as a ID-resonant tunneling device. The reflec- 
tion and transmission amplitude through the dot can be written in terms of the scattering 



17 



I 

o 



4> 



4.6 - 



4.4 



4.2 - 



4.0 




■5 5 10 

Magnetic Field, B [mT] 

FIG. 11. Oscillatory part of the ring resistance as a function of the magnetic field. Measured 
for a two-terminal AB ring with an artificial impurity embedded in one of its arms. Taken from 



Ref. 118 



amplitudes of the adjacent tunneling barriers and a phase 9 accumulated by the motion in 
the dot. A variation of the plunger voltage at the dot is simulated by a variation in 9. 

The predicted amplitudes of the h/e and the h/2e harmonics are shown in Fig. [10]. The 
h/e component grows as a conductance peak is approached from the left. Near the peak 
the h/e component drops and vanishes for some value of 9. The current oscillations are 
then dominated by the h/2e oscillations. Upon further increase of 9, the h/e component 
changes sign and grows again. The change in the sign of the h/e component is detected 
experimentally as a jump by 7r in the AB phase. 



Jacoby et al. Jl8| corroborated their theoretical arguments by a measurement of the AB 
oscillations through a ring with an artificial impurity. A small gate inserted in one arm of 
the ring creates a potential barrier for the electron motion. The potential at the gate can 
be varied by means of a metallic air bridge. The impurity replaces the quantum dot in the 
earlier experiment || . The AB oscillations of the ring resistance are displayed in Fig. |TT] for 
increasing impurity strength. One observes that the h/e component has a rigid phase. The 
h/2e oscillations become strong when the h/e oscillations vanish. The phase jump in the 
h/e oscillations is thus associated with a qualitative change in the AB oscillations: The h/e 
oscillations vanish and the AB oscillations show h/2e periodicity. 

The symmetry ( |3.1|) of the two-terminal conductance strictly holds only in the linear 
response limit of infinitesimally small currents and voltages. Bruder et al. [I?]] investigated 
the nonlinear response regime using a tunneling Hamiltonian description of the quantum 
dot. They expressed the current through the AB ring in terms of Green functions of the 
quantum dot. This description includes charging effects in a non-perturbative way. Bruder 
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et al. recover a symmetric conductance in the linear response limit. In the nonlinear response 
regime they find derivations from the symmetry under sign change of the external flux. As 
a consequence, the phase of the AB oscillations can change continuously with the voltage at 
the quantum dot when a finite voltage difference is applied across the AB ring. 



B. Aharonov Bohm current 

The phase of the AB oscillations in a two-terminal measurement is rigid and yields lit- 
tle information about the transmission phase through the quantum dot. The latter phase, 
however, strongly affects the amplitude of the AB oscillations as was first shown by Hack- 
enbroich and Weidenmiiller ||16| , [19| . Their approach starts from a description of the dot and 
the AB ring in terms of a tunneling Hamiltonian. With little modification this Hamiltonian 



also serves as the starting point for the discussion in Sec. [TV} We introduce the model in 



Sec. [I I IB 1| . The derivation of the AB amplitude is sketched in Sees. 1 1 1 B 2| - |HI B~3[ The 



calculation is based on single-particle scattering theory and the Landauer formula. This re- 
stricts the approach to the regime near the conductance peaks and to temperatures kT -C A, 
where A is the single-particle level spacing of the quantum dot. 



1. Tunneling Hamiltonian 



Consider the system schematically represented in Fig. |12| A conducting ring threaded by 
the magnetic flux <3> is connected to two leads. The ring and the leads may support several 
transverse channels. A quantum dot is embedded in one arm of the ring. The quantum dot 
is weakly coupled to the ring by potential barriers on either side of the dot. The quantum 
dot is in the Coulomb blockade regime. The electrostatic potential at the dot is controlled 
by the gate voltage V g . 

The system is described in terms of the tunneling Hamiltonian 

H = H + H T , (3.6) 

where H describes the isolated subsystems and H? the couplings between these subsystems. 
Explicitly, H is given by 

H = £ / dEEa^ E a r mE + £ e^d, + £ E x c\c x + U, (3.7) 

mr J i A 

where r = 1,2 labels the two leads, m the channels in either lead and i and A the single- 
particle states in the AB ring and the quantum dot, respectively. The respective single- 
particle energies are denoted by E (the longitudinal energy in a channel), and £\. The 
interaction U has the form 
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L V R 

FIG. 12. An AB ring threaded by the magnetic flux $ is connected to two external leads. A 
quantum dot is embedded in one arm of the ring. Tunnel barriers between the dot and the ring 
are modeled by matrix elements Vl,r- The coupling between the ring and the leads is described 
by the matrix elements W\ 2 . 



U 



^o(iV 2 



N), 



(3.8) 



where N = J2x c x c x * s the total number of electrons on the dot, Uq = e 2 /C the charging 
energy, and C the total capacitance between the dot and its surroundings. We use the 
standard picture of the Coulomb blockade |23|] and assume that the energies £\ vary linearly 
with the gate voltage 



c x 



9- 



(3.9) 



where a is some function of the capacitance matrix of the system. Below we assume that 
the energies E\ are non-degenerate. 
The couplings have the form 



H T = (£ J dEWUEW^d, + H.c.) + £ Kcfa + H.c. 

mri ipX 



(3.10) 



where the matrix elements W couple states in the ring to states of the leads, V provides the 
coupling between states in the ring and the dot, and p = L,R labels either side of the dot. 
All matrix elements do not change appreciably on the scale of the charging energy. This 
weak energy dependence will be neglected below. 

It follows from time-reversal invariance that for vanishing magnetic flux through the 
ring, all matrix elements can be taken real. To account for non-zero flux $ through the AB 
ring, we attach to each matrix element V^- a factor exp(z</>), where = 27r$/$ - We can 
thus write 



Xi 



V£exp(-i<j 



V x fexp(i<p) = v 



R 
Xi- 



(3.11) 
(3.12) 
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where v Xi and v Xi are real. This parameterization is adequate whenever the single-particle 
states both in the AB ring and in the dot do not change appreciably with flux, i.e. whenever 
the flux through each arm of the ring and through the dot is smaller than $ - 



2. Scattering matrix 

Near transmission resonances and for temperatures kT much less than the single-particle 
level spacing, scattering through a quantum dot may be described in a single-particle picture. 
The scattering matrix may then be derived explicitely. We first sketch the caclulation for 
noninteracting electrons (the limit of zero charging energy Uq — ■> 0) and then indicated the 
changes resulting from nonzero Uq. 

The scattering matrix 5 = 1 — 1-niT at energy E and (dimensionless) flux <fi is obtained 
from the Lippmann-Schwinger equation 

T = H T + H T (E - H + z?7) _1 T (3.13) 

for the transition operator T. Here, rj is positive infinitesimal. Iteration of the Lippmann- 
Schwinger equation yields the Born series 

T = H T + H T — H T + H T — ———H T — — ■ —Hx H . (3.14) 

E — Hq + 17] ti — H + 17] h — Hq + 17] 

For non-interacting electrons the right-hand-side reduces to a geometric series which can 
easily be resumed. For the S-matrix element S™ n connecting channel n in lead s with 
channel m in lead r one obtains the result 

S™ n (E, 0) = 5 rs 5 mn - 2m[WD;^]Z n - (3.15) 

where we used the shorthand W and V = V L + V R for the coupling matrix elements and 
introduced the propagators 

{E - ei)5 lk + i7r[W^W} ik , (3.16) 
D dotM — (E — S x )5 Xfl - [VD^V^, (3.17) 

for the ring and the dot, respectively. The partial width amplitudes 7^ A , 7^ A are given by 

7 ^ A = V2^[W(D lins )-W} r mX , (3.18) 
l s n , = ^[W{D\ mg )- l V^. (3.19) 



The terms on the right-hand-side of Eq. (|3.15| ) have a simple physical interpretation. The 



first term involving the Kronecker deltas describes the reflection from channel m back into 
the same channel. The second term is flux and voltage independent and accounts for the 
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scattering of elctrons through the free arm of the ring. The last term describes the scattering 
of electrons through the dot. The entire flux and voltage dependence resides in this term 
(via both D dot and the partial width amplitudes). Note that the flux dependence does not 
reduce to a phase factor exp(i0) as one naively might have expected. The origin for the 
more complicated flux dependence is the multiple scattering of electrons through the ring 
which gives rise to higher harmonics of the current oscillations. 

In the Coulomb blockade regime, the quantum dot is weakly coupled to the AB ring. 
The quantum dot resonances are then isolated i.e. their spacing is much larger than their 
width. In this case and for E close to a resonance, one can approximate -Ddot by the diagonal 
matrix 

(-Ddot) A/i = (E — £\ — A£ x + zr A /2)V (3.20) 

Here, AS\ is the energy shift (of the resonance position with respect to the bound state 
energy of the isolated system), and T\ the total width of the resonance with index A. Using 
the unitarity of S, one can show that for isolated resonances 

Ta = EI7LJ 2 = EI7™a| 2 - (3-21) 

mt mt 

Note that because of the oscillatory behavior of the partial width amplitudes, also T\ is an 
oscillating function of flux. We further note that the S'-matrix ( |3.15| ) has the well-known 
symmetry properties 

5(0)5t(0) = &(<j>)S(<j>) = i, (3.22) 
S(0)S*(-</>) = Sr(<f>)S{-4)= 1, (3.23) 

that follow from current conservation and full time-reversal invariance (including a reversal 
of the magnetic field). Both symmetries together imply that the linear conductance is an 
even function of magnetic flux plj and impose the phase rigidity as discussed in Sec. [Ill A. 



The S-matrix ([3.15 ) for non-interacting electrons may be generalized to the case of 



nonzero charging energy using the Hartree-Fock approximation. For isolated quantum dot 



resonances, one finds that the scattering matrix S(Ep, 0) at the Fermi energy has the 
matrix elements 



]rs 
mn 



S™ n (E F , <f>) = S rs 6 mn - 2m[WD^ g W^}l 

-i^2Y m x(EF -E x + zr A /2)- 1 7ll, (3.24) 



where 7 m A, 7nA, and T x are defined in Eqs. (|3.18|) , ( [3.19D , and (|3.21|) , respectively. The S- 



matrix has the same form as in the non-interacting case ( |3.15| ). However, the single-particle 
energies E\ have been replaced by the resonance energies E\ which are determined by the 
set of self-consistent equations 
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E x = S x + AS\ + f/oE M> (3-25) 
where the the average occupation probability (n\) of the level A is given by 

(n x ) = - [°° dEf(E - E F ) /9 , 2 - (3-26) 

n J-oo \E — E\ + lY \/ Z\ A 

Here f(E—Ep) = [l+exp((£'— i?i?)/fcT)] _1 denotes the Fermi function. The solution to these 
equations is easiliy found when the single-particle levels are non-degenerate. Whenever one 
level becomes filled, all higher levels are shifted upwards by an amount equal to the charging 
energy. In effect, the Hartree-Fock approach reduces to a picture of isolated resonances with 
a 'stretched' level spacing due to the Coulomb interaction. 

We note that the mean-field approach used here is only valid close to the resonances 
and at temperatures kT much less that the single-particle level spacing A in the quantum 
dot. The conductance at higher temperatures kT ~ A can no longer be reduced to a single- 
particle scattering problem. The sequential (incoherent) tunneling through a quantum dot 
in the regime kT ~ A has been studied by Beenakker p3j and by Meir et al. [ j34j] . 



3. Conductance 



The dimensionless conductance g = (h/e 2 )G of the AB ring with the quantum dot is 
obtained from the multi-channel Landauer formula 



g = 2 dE l-^L) V \t mn (E)\ 2 . (3.27) 




Here t mn (E) = S^ n (E) is the transmission amplitude through the ring for an electron 
entering the ring via channel n in lead two, and leaving it via channel m in lead one. The 
derivative of the Fermi function / is given by -(df/dE) = (AkT)- 1 cosh -2 ((E - E F )/2kT), 
and Ep is the Fermi energy in the leads. A factor two accounts for the spin degeneracy of 
the electron. 

It is assumed that both the charging energy Uq and kT are much larger than the resonance 
widths T^. Moreover, kT shall be much smaller than the single-particle level spacing A. 
Then an appreciable current can pass the dot only if a resonance in the dot is close to the 
Fermi energy, Ex ~ Ep- The contribution of other resonances to the transmission amplitude 
can be neglected. To simplify notation, we suppress the index A below. According to 
Eq. ( |3.24[ ), the transmission amplitude has the form 



y y* 

L mn — (ring, mn E — E' -\- iY / 2^ \<->.£0) 
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where E' is the resonance energy and t r in g ,mn = —^i^J2ik ^mi(^°)ik^nk * s ^ ne transmission 
through the free arm of the ring. This contribution is independent of <fi and V g and only 
weakly dependent on energy (this energy dependence will be neglected). 

We now confine ourselves to a symmetric dot where v\ = vfexpiix) = v%- Due to time- 
reversal symmetry, Vi can be chosen real, and x can only take the values and it. Substituting 
Eq. Q3.28| ) into the Landauer formula, summing over the channels and integrating over 



energy, one finds that the conductance can be written as 

f 

9 = 0ring + [1 + COS(0 - X)]X — 5 SUl(£ - f3) 

+ [1 + cos(0 - x)]y^~ cosh- 2 {^^\ , (3-29) 

where g T i ng , x, y are positive coefficients, £ is a real phase shift, and T = (1 + cos[0 — x]) _1 r 
is the flux-independent width. We note that g is an even function of the magnetic flux 
<p as required for a two-terminal measurement. Both the amplitude B and the phase (3 
are functions of (E F — E')/{kT) (see Fig. The phase (3 may be identified with the 

transmission phase through the quantum dot. It takes the value 7r/2 at E = E', while 
approaching for E — > — oo and 7r for E — > oo, respectively. 

Equation ( |3.29| ) allows us to discuss the current through the AB ring as a function 
of energy (voltage V g ) and magnetic flux. We refer to the first term on the right hand 
side of Eq. (|3.29|) as background term, to the second term as interference term, and to 
the third term as resonance term. The background term comprises the bulk part of the 
current. This term is independent of both energy and flux. It is due to transmission of 
electrons through the free arm of the AB ring. This term is subtracted in the experimental 
analysis from the total current through the device. The other terms are smaller than the 
background term by a factor T/kT <C 1. The interference term is due to the interference of 
the background amplitude with the transmission amplitude for passage through the quantum 
dot. This term displays resonant behavior and depends explicitly on the geometry of the ring 
through the phase shift £. The amplitude of the resonance term has the usual temperature 
dependence known for thermally broadened resonances [^| of quantum dots directly coupled 
to leads. However, in the AB geometry this term is flux-modulated and contributes to the 
AB oscillations. 

We now focus on the AB phase. As a function of the voltage at the dot, only the 
amplitude of the current oscillations is changed. The AB phase is unaffected unless the 
interference term changes sign. In this case, the AB phase jumps by tt. The location of the 
phase jump depends on the system specific phase shift £. Eq. (|3.29|) predicts that all higher 
harmonics vanish identically. This is a consequence of keeping only the resonance closest to 
Ep and neglecting the far away resonances. Such resonances as well as an asymmetric dot 
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FIG. 13. Amplitude of the second term of Eq. ( |3.29D as a function of {Ep — E')/kT for 
£i = 7r/2. The dashed line shows the phase extracted from the sign of i3sin(£i — /?). Taken from 



Ref. p|. 



produce higher (e.g. $o/2 periodic) harmonics which may be observable whenever the $o 
periodic harmonic vanishes. 

Inelastic scattering in the quantum dot destroys the unitarity of the 5-matrix. Phe- 
nomenologically, this increases the widths of the transmission resonances T by the addition 



of an inelastic width p5,pT and thereby reduces the amplitudes of both the interference 
term and the resonance term. The fundamental property of the current to be an even func- 
tion of the applied flux remains unchanged. As a result, with properly rescaled amplitudes, 
Eq. ( |3.29|) accounts for the current through the system even in the presence of inelastic 
scattering. The basic conclusions concerning the temperature and flux dependence of the 
various terms remain valid. 
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IV. TRANSMISSION PHASE THROUGH A QUANTUM DOT 



The experiment of Yacoby et al. |J proved for the first time that part of the transmission 
through a quantum dot in the Coulomb blockade regime is coherent. The phase of the 
transmission amplitude through the dot could not be measured in this experiment due to 
the phase rigidity discussed in Sec. |TJ. The problem of phase rigidity was solved by Schuster 
et al. [0] who replaced the two-terminal AB ring by a multiple-terminal device. The AB 
phase in the Schuster-experiment increased roughly by tt whenever the gate voltage at the 
quantum dot was swept through a transmission resonance, and the profile of the phase 
evolution was well described by a Breit-Wigner formula. At the same time, the phase 
displayed unexpected properties: (i) It was nearly identical for all resonances investigated in 
the experiment, and (ii) between each pair of adjacent resonances the phase sliped roughly 
by tt on a very small energy scale. Both observations (i) and (ii) were totally unexpected. 
The theoretical work addressing these observations is reviewed in this section. 

Our understanding of the transmission phase is largely based on quantum mechanical 
scattering of non-interacting particles. The one-dimensional analogue of a quantum dot is 



a double-barrier well as depicted in Fig. |l4](a). The tunneling barriers mimic negatively 
biased gates on either side of the dot. A change in the gate voltage at the quantum dot is 
modeled by a shift of the potential in the well. The calculation of the transmission amplitude 
t through the double-barrier well is a textbook problem. The transmission coefficient \t\ 2 
and the transmission phase arg(i) are shown in Fig. 0(b) as a function of the well potential. 
The phase increases smoothly by tt over the width of a transmission resonance and stays 
nearly constant between resonances. As a result, the phase of neighboring resonances always 
differs by tt. We refer to such resonances as off-phase resonances below; resonances with 
similar transmission phase will be called in-phase resonances. The appearance of off-phase 
resonances is not a special property of the double-barrier quantum well but holds generally 
for scattering through time-reversal invariant strictly one-dimensional systems. In such sys- 
tems, the wave functions can be chosen real. The wave functions representing neighboring 
resonant states differ by one additional node, and with each node the phase of the trans- 
mission amplitude increases by tt. As a result, the phase evolution in the one-dimensional 
model is qualitatively different from the phase evolution measured in Refs. 0,0). Thus two 
questions arise: First, how can in-phase resonances occur in the Schuster- (and Yacoby-) 
experiment, and why are all resonances in phase? Second, why does the phase slip between 
resonances and why so sharply? 

Although we will focus on the transmission phase, we note that the underlying questions 
go much deeper and concern the nature of electron transport through quantum dots. Starting 
with the work of Jalabert, Stone and Alhassid in the early 1990's, the theoretical work 
on quantum dots has almost exclusively addressed universal aspects described by random 



26 




FIG. 14. (a): Double-barrier well potential, (b): Transmission coefficient |t| 2 (top) and trans- 
mission phase arg(t) (bottom) both as a function of W. Parameters are Vq = 10.0, E = 2.0. All 
energies are given in units of fo 2 /2ma 2 . 

matrix theory. Employing random matrix theory, quantities like the statistical distribution 



of the conductance peak heights |26[ and parametric peak-height correlations p7|>eH nave 
been calculated and tested in various experiments. Because of the enormous success of these 
calculations random matrix theory has emerged as a paradigm for the theoretical description 
of quantum dots. However, in recent years a number of experiments showed clear deviations 
from random matrix behavior. Among these are transport experiments on small vertical 



quantum dots [29|, conductance peak-height measurements on lateral quantum dots [30 
and the Coulomb-blockade interference experiments |j],[7|. We will show below that the 
failure of random matrix theory for the experiments |||7| may be traced back to features of 
the short-time dynamics in the quantum dot. Under certain conditions such non-universal 
features can totally dominate the quantum transport in the tunneling regime and generate 
strong deviations from random matrix predictions. 

The material covered in this section is organized systematically so that papers with re- 
lated theoretical ideas are discussed together. This ordering displays the underlying physical 
concepts most clearly. Sometimes the systematic ordering could only be achieved at the ex- 



pense of giving up a strictly chronological order. In Sec. IV A we summarize the scattering 
theory through a region of interacting electrons, and derive the transmission phase through 
a weakly coupled quantum dot. This section forms the basis for the arguments presented in 



the subsequent sections. In Sec. |IV B| , we discuss a general theorem known as the Friedel 



sum rule. The Friedel sum rule relates the amount of charge added to a conducting system 
to the sum of all scattering phase shifts. We clarify the relation between the scattering phase 
shifts that enter in the Friedel sum rule and the transmission phase measured in a quan- 
tum dot interference device. The phase lapse between resonances is discussed in Sec. [IV C 
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Mechanisms for in-phase resonances are summarized in Sec. [IV D| . We emphasize, that de- 
spite of the arguments for in-phase resonances, non of the theoretical studies excludes the 
possibility of out-off-phase resonances. Further experimental work is necessary to clarify 
whether the latter may or may not be found in quantum dots. 



A. Transport theory 



is 



The mean-field theory for transport through a quantum dot described in Sec. [Il l B 
only valid near the conductance peaks and at temperatures much smaller than the mean 
level spacing A. When these conditions are not met, transport through a quantum dot can 
no longer be reduced to a mean-field problem. Several techniques have been developed to 
calculate the current through interacting electron systems. Our presentation follows the 



approach of Refs. |31~]-|3~3|| which is based on earlier work of Meir, Wingreen, and Lee fl34| , |2"4" 



1. Generalized Landauer formula 

We consider an Aharonov-Bohm ring with a quantum dot embedded in one arm. The 
arm containing the dot is modeled by the tunneling Hamiltonian 

H = H L + H R + H QD + H T , (4.1) 

H L{R) - 2^ € k a k a k > l 4 ^J 

k 

Hqb = E &4c A + \u (N 2 - N) , (4.3) 

Ht = Y: (v k L x a L k c{ + H.c.) + L <- R . (4.4) 

k,\ 

Here, Hi and Hr, respectively, describe the regions to the left and right of the QD, Hqd is 
the Hamiltonian of the isolated QD including the charging energy, and Ht represents the 
tunneling of electrons in and out of the QD. All energies are counted from the Fermi level 
in the leads. For simplicity, we assumed only one transverse channel in the left and right 
region. In the absence of a magnetic field, the matrix elements V k ' x can be chosen real. 

An exact formula for the current through the ring can be derived using the nonequilibrium 
Keldysh formalism [pi . The current is expressed in terms of the Fermi function in the leads 



and local properties of the interacting region. In general, the formula for the current includes 
inelastic scattering, spin flips, and even scattering processes of several electrons. However, 
in the linear response limit and for sufficiently low temperature only elastic processes are 
allowed by energy conservation. Then, the conductance through the Aharonov-Bohm ring 
reduces to the generalized Landauer-type formula G ring = (2e 2 //z)T ring , where the total 
transmission probability T ring is given by 
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t(£)exp[27ri$/$ ]| 2 . (4.5) 

Here to and t(E) denote the transmission amplitude through the free arm and the arm 
with the quantum dot, respectively, $ is the magnetic flux through the ring, and $o — h/e 
the elementary flux quantum. Note that we have neglected higher harmonics. From the 
interference term in Eq. (|4.5| ) one can extract the amplitude of transmission through the 
arm with the quantum dot 

t QD = J dE(-df/dE)t(E). (4.6) 

We identify #qd = arg(tQD) with the transmission phase through the quantum dot. The 
transmission coefficient through the arm with the dot is given by 

Tqd = / dE(-df/dE)\t(E)\ 2 . (4.7) 

Both 6qd and Tqd can be measured in quantum dot interference experiments. The trans- 
mission amplitude t(E) can be expressed in terms of the retarded Green function G Xfl of the 
dot 

t(E) = £ V x L (E)G r x ^E)V R *(E) , (4.8) 
\,fi 

where we introduced the tunneling amplitudes V X L,R (E) = \2np L ' R {E)] l / 2 V^ R - ) x with the 
density of states p L ' R (E) in lead L, R. To simplify notation, we will drop the argument E 

L Ft 

and write V x ' below. We emphasize that the Green function G r Xfl must be calculated in 
the presence of interactions and tunneling. 

In the limit to = of vanishing transmission through the free arm, the ring conductance 
reduces to the two-terminal conductance Gqd = (2e 2 /h)Tqp through a quantum dot coupled 
directly to reservoirs. We derive Gqd for the case when the Green function is diagonal in 
the single-particle basis of the dot. Substitution of Eq. (4.8) into Eq. (|4.5|) then yields a 



product of two Green functions which can be evaluated using the relation 

1 Tm C r 

G r xx Gl* x = — (G xx - G r xx ) = — (4.9) 

^aa — ^aa J-rnz^AA 

where S is the self-energy. The diagonal elements Saa = ^(T^ + Y x )/2 are given in terms 
of the partial widths T X ,R = \V X ' R \ 2 for decay of the state A into the left and right lead, 
respectively. The current takes the form 



This result was first obtained in Refs. [34,p4"|. Note that —(l/ir)lmG xx (E) is the local 



level density in the presence of interactions and tunneling. The above derivation can be 
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generalized to include the electron spin. In this case the conductance takes a form similar to 
Eq. (14.101) , but the factor 2 is replaced by an explicit summation over the electron spin. The 
resulting expression has been used to study the conductance in the Kondo regime PHE?3 . 



2. Retarded Green function 



To derive the dot Green function we employ the equations-of-motion method pi |. The 
derivation starts from the retarded Green function G r Xfl (t), defined by 

Gyt) = -i0(t)<{ CA (t),<(O)}}, (4.11) 

where the curly brackets denote the anticommutator. The expectation value is the ther- 
mal average with respect to the Hamiltonian H. The operator c\(t) is the solution of the 
Heisenberg equation 

i^ = [c»,fl]. (4-12) 



Differentiating Eq. fl4.11|) with respect to t and substituting Eq. (|4.12|) one obtains the 
equation of motion 



d_ 

Of 



G r Xfl (t) = -i5(t)8 x ,-9(t)({[c x ,H],cl}). (4.13) 



When H is quadratic in the particle operators, one can express the right-hand-side of this 
equation in terms of two-particle Green functions. This yields a closed set of equations which 
is readily solved for the exact Green functions. Since the quantum dot Hamiltonian H is 
not quadratic, one generates higher-order Green functions which must be approximated to 
obtain a closed set of equations for G r x (t). From the solution one finds the Green function 
G r Xfl (E) upon Fourier transformation. It is illustrative to study G r x ^(E) in various limiting 
cases. 

L Ft 

Isolated dot: This case is obtained for zero coupling V x ' =0. Since no electrons can 
tunnel into or out of the quantum dot, the total number of dot electrons N is conserved. The 
Fock states of the dot can be separated into classes with N electrons. The grand-canonical 
expectation value of an operator A can then be expressed in terms of canonical expectation 
values 

oo 

(A) = J2 Pn(A) n , (4.14) 

7V=0 

where the equilibrium probability Pjy to find N electrons on the dot is given by 

tr N exp(—(3H) 



JV 



E^ =0 tr M exp(-/3iO 



(4.15) 
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The retarded Green function follows upon substitution of A = —i$(t)c\(t)d (0) and A 
— i9(t)c' (0)cx(t). Fourier transformation yields the result 



N=0 



1 " (nx) 



N 



E-(£ x + N- U ) 

(hx)N 



iS 



E - (E x + (N - 1) • U + i5 



(4.16) 



where 5 — > + . Note that each level A contributes twice: The first term on the right-hand- 
side of Eq. ( 4.16 ) accounts for electron propagation through the level A while the second term 
describes the propagation of a hole. The particle-hole structure originates from the charging 
energy which generally assigns the addition or removal of electrons a different energy. We 
note that the particle— hole structure disappears for Uq — > 0. 

Noninteracting electrons: The limit Uq — > corresponds to noninteracting electrons. 
In this case the equations of motion close and one readily finds 
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(4.17) 



The self-energy T,(E) = Y, L (E) + Y, R (E) has two contributions resulting from the coupling 
to the left and right lead. They are given by 
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(4.18) 



Since the states in the leads are dense, one can replace the sum in Eq. ( |4.18j ) by an integral. 

L Ft 

The imaginary part of this integral yields the widths T x ' for decay of the state A into the 
left and right lead, while the real part gives an energy shift. For weak coupling r x ' R A 
and close to the resonances, the Green function can be approximated by the diagonal terms 
G x^ = yielding 



GUE) = 5. 



E-Sx + irx/2' 



(4.19) 



with the total width Tx — ± A 
single-particle scattering theory. 



In Sec. Ill B we arrived at the same result using 



3. Weak-coupling limit 



The weak-coupling regime is characterized by Tx kT, A. In this regime the Green 
function can be approximated by 
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(4.20) 



Note that due to the coupling all quantum dot states acquire a finite width. The widths are 
identical to the widths obtained in the non-interacting case. Corrections arise for T\ > kT 
and lead to a normalization of the self-energy. In particular, the approximation ( |4.20| ) fails 
in the Kondo regime where higher iterations of the equations of motion become important. 

In the weak-coupling limit both the transmission amplitude and the transmission co- 
efficient can be obtained analytically. We first substitute the weak-coupling result for the 
Green function (|4.20|) into Eq. ( f4.8|) . The transmission amplitude (|4.6| ) and the transmis- 
sion coefficient ( [4.7|) then reduce to sums of integrals which can be computed using contour 
integration. For the transmission amplitude, the integrals are of the type 
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E-S + iT/2 2-ni 
and for the transmission coefficient they have the form 
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Here = 1/kT is the inverse temperature, / = 1/(1 + exp((3E)) is the Fermi function, and 
ip^(z) the trigamma function | 36| . Combining results, the transmission amplitude reduces 
to a sum over energy levels and occupation numbers, 
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(4.23) 



The transmission coefficient Tqd and the conductance Gqd = (2e 2 //i)TQD reduce to similar 
sum. The result takes a simple form in two limiting cases. First, in the cotunneling regime 
deep in the transmission valley iV we have Pm — $m,n- The first term in Eq. ( |4.22j ) vanishes 
exponentially, while the trigamma function can be approximated by its asymptotic expansion 
The conductance then simplifies to 
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(4.24) 
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This is the standard result |37j for the conductance in the cotunneling regime. Second, 
near the transmission resonances, we can neglect the terms involving the trigamma function 
which are smaller than the Fermi function by a factor /3T\. After a little algebra one finds 
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(4.25) 



This is the well-known result for the conductance in the regime of sequential tunneling that 



was obtained independently by Beenakker p3[ using a master equation approach and by 
Meir, Wingreen, and Lee |53J using the generalized Landauer formula. 

The results derived in this subsection solve the problem of electron transmission through 
a weakly coupled dot. They express the transmission phase in terms of the coupling matrix 
elements and the single-particle energies of the quantum dot. These quantities must be 
obtained from a dynamical model of the dot and its coupling to the leads. Various such 
models are studied below. Deviations from the weak-coupling regime are addressed in 
Sees. [TVTJ^, [TyTJg. 



B. Priedel sum rule 



Levy Yeyati and Buttiker |0| were the first authors to study the origin of the in-phase 
resonances. They discussed the subject in terms of the Friedel sum rule [^,^]. The Friedel 
sum rule is a powerful relation that links the number of localized electron states in a con- 
ducting region to the sum of all scattering phase shifts at the Fermi energy. The sum rule 
holds even in the presence of interactions |38| . There is an analogue of the Friedel sum rule in 
scattering theory. Known as Levinson's theorem it relates the number of bound states 
in a scattering region to the sum of the scattering phase shifts at zero energy. Employing 
the Friedel sum rule, Levy Yeyati and Buttiker argued that the addition of a charge AQ 
to a conducting region Q should change the phase of the transmission amplitude t through 
that region according to 



AQ/e = Aarg(t)/vr. 



(4.26) 



If the region Q is chosen to include only the quantum dot, then the addition of one electron 
to the dot should cause a phase increase by n. Consequently, successive resonances would be 
off-phase by tc rather that in-phase as in the experiments |||7|. Levy Yeyati and Buttiker 



emphasized, however, that the charge AQ entering in Eq. ( f4.26|) is not the additional charge 
on the quantum dot but rather the additional charge within the whole coherence volume. 
This volume includes the quantum dot and the AB ring around the quantum dot. A variation 
of the gate voltage at the quantum dot may also modify the electrostatic potential in the 
vicinity of the dot. This could result in the addition of extra charge to the AB ring outside the 
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FIG. 15. (a) In strictly ID systems, even and odd resonant levels alternate in energy, (b) In 
quasi-lD or higher dimensional systems, they do not necessarily alternate, allowing for sequences 
of resonances with the same parity. 

quantum dot. Levy Yeyati and Biittiker studied a situation where with each electron added 
to the quantum dot, the extra charge 5e is added to the AB ring. In a model calculation 
for a two-terminal device using the values 5 ~ 0.3 and S ~ 1, sequences of two or three 
consecutive in-phase resonances were found. 

The connection between the Friedel sum rule and in-phase resonances was reconsidered 
and clarified by Lee [El] and Taniguchi and Biittiker |42| 



The Friedel sum rule makes a 
statement about the sum of all scattering phases. The sum rule can be stated in the form 



AQ/e = [AlnDet(S)]/2vri, 



(4.27) 



where S is the scattering matrix. Writing the eigenvalues of S in the form e 2m % the right 
hand side represents the change in the sum ctj-n = J2i=i ctij-n of all eigenphases, where N is 
the dimension of S. Thus, 



AQ/e = Aa/n. 



(4.28) 



The transmission amplitude t, on the other hand, is a matrix element of S in a basis of left 
and right moving scattering states. In general, the phase of t is not simply related to the 
sum of the eigenphases Aa ^ Aarg(t). The relation (|4.26| ) is thus different from the Friedel 
sum rule ( |4.27| ), ( |4.28[ ). Generally, the expression ( |4.26| ) is not expected to be correct. 

The difference between Eqs. ( |4.26|) and (|4 . 2 8|) can be explicitely demonstrated in the 
scattering through a quasi-lD system |fj"l| . We assume a mirror reflection symmetry x — > 
—x and the absence of a magnetic field. The scattering region at |x| < R is connected 
to two single-channel leads at x < —R and x > R, respectively. The scattering states 
can be decomposed into even and odd scattering states. For |x| > R these states read 
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Note that there are 
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two scattering phases a e , a a corresponding to two scattering channels, one in each lead. 
Rather than in a basis of even and odd scattering states, one can formulate the scattering 
problem in a basis of right and left moving scattering states, ipi(x) = [if) e (x) — ip (x)}/2, 
i/j (x) = [i(j e (x) + ip (x)}/2. The scattering matrix in the new basis takes the form 




S^ ir ,j, (4.29) 
where t, t' are the transmission amplitudes and r, r' the reflection amplitudes from the left 



and right, respectively. From the relations between the basis vectors and Eq. ( 4.29 ), one 
obtains the transmission amplitudes 

t = tf = ie ia sin (3 (4.30) 

in terms of the phases a = a e + a Q and (3 = a e — a . In true ID systems, even and odd 
resonant states alternate in energy (see Fig. p!5](a)) and the angle (3 is restricted to the range 
< (3 < 7r. Then both the total phase a and the transmission phase arg(t) change by 7r 
between neighboring resonances. In quasi- ID systems, however, even and odd states do not 
necessarily alternate in energy, and there may be sequences of states with the same parity 
(see Fig. 0(b)). Within such sequences, (3 varies by more than n and sin (3 changes sign from 
one resonance to the next. Neighboring resonances then have the same transmission phase 
arg(t) while the total phase a changes by 7r in keeping with the Friedel sum rule Q4.27|) , 
( |4.28| ). Note that Eq. (|4.30| ) predicts a transmission zero (sin/? = 0) between resonances 
with the same parity. A transmission zero corresponds to a singular point of the transmission 
phase. The phase jumps abruptly by it when the system is swept through the transmission 
zero. 

The arguments presented above demonstrate that in-phase resonances do not violate 
the Friedel sum rule. They also show that the spatial dimension of the scattering region is 
important: While neighboring resonances are always out-off-phase in strictly ID systems, 
both out-off-phase and in-phase resonances can be found in higher dimensional systems. 



C. Phase lapse 

1. Breit-Wigner model 

The sharp phase lapse between resonances was investigated by Hackenbroich and Wei- 
denmiiller |4^| for the weak-coupling regime V < kT < A. Similar to the calculation for the 
two-terminal case, cf. Sec. [Ill B| , the multi-terminal transmission amplitude was obtained 
using scattering theory and the diagonal approximation for the dot propagator. In this 
approximation the transmission amplitude reduces to a sum of Breit-Wigner resonances. 
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The internal dynamics of the dot and the coupling to the leads enter via the resonance 
energies E\ and the partial width amplitudes V^ m for the decay of the resonance A into 
channel m of the lead p. Hackenbroich and Weidenmiiller did not calculate the resonance 
parameters from a dynamical model. They rather assumed these parameters to be identical 
for all resonances and investigated the consequences of this assumption for the transmission 
phase: By assumption, all resonances showed in-phase behavior. At each resonance the 
phase increased by 7r on a scale ~ kT. Between the resonances, the phase sliped by 7r on 
the scale of the intrinsic resonance width T. The Breit-Wigner model (|3| thus predicts a 
phase lapse between resonance peaks on the scale of the intrinsic resonance width T. 

We note that the experiment [0] was not performed in the weak-coupling regime. Rather, 
the intrinsic resonance width T in the experiment was estimated to be of the order T ~ 
(3 — 4)kT. The Breit-Wigner model of Hackenbroich and Weidenmiiller therefore cannot 
explain the phase lapse observed in 0. Moreover, there is theoretical evidence (see Sec. 
|I V C 2|) that the diagonal approximation for the quantum dot propagator fails between the 
resonances. As a result, the transmission phase may display an sharp lapse rather than a 
smooth decrease on the scale T as predicted by the Breit-Wigner model. 



2. Transmission zero 



In 1998, a number of authors |33-fT7| studied models for the transmission through a quan- 
tum dot and found abrupt jumps of the transmission phase. The jumps occur ed at singular 
points where both the real and the imaginary part of the complex transmission amplitude 
vanished. Transmission zeros were found in numerical studies exploring Fano resonances 



44|-|46[ and the interplay of multiple resonances |47j. In Refs. |44j-[46| the quantum dot was 
modeled as a quasi-lD or 2D region of regular shape. In all studies time-reversal invariance 
was assumed, the dot was connected to two single-channel leads and electron-electron inter- 
actions were neglected. In-phase resonances occurred either due to sequences of states with 
the same parity or by assuming J|7j] state-independent coupling matrix elements to 



the leads. A more general approach to the connection between transmission zeros and in- 
phase resonances was presented by Lee (see the discussion after Eq. (|4.30|) ). He showed 
that the transmission always vanishes between neighboring in-phase resonances, provided 
(i) the scattering region is connected to two single-channel leads and (ii) the system is 
time-reversal invariant. 

Why do the models [^-|47|| display an abrupt phase jump and not a smooth phase slip 
on the scale T as predicted by the Breit-Wigner model [53|? This question was addressed 



by Sun and Lin [[17]] for a quantum dot with symmetric tunneling barriers. It turned out 



that the diagonal approximation for the quantum dot propagator used in Ref. [43] fails 



for time-reversal invariant systems coupled to single-channel leads. Here we present an 



36 



argument that proves this failure also for asymmetric barriers. Consider a quantum dot 
that is coupled to two single-channel leads, denoted L, R. We assume fcT-cA and include 
electron-electron interactions in a mean field way. The scattering matrix of the system can 
be derived with the methods of Sec. |111 BL For the scattering amplitude S mn at energy E 



one has 

S mn (E) = 5 mn - 2inJ2W^(D~ 1 (E)) tMU W nU} (4.31) 



where m, n refer to the physical channels. The indices /i and v of the inverse propagator 
D^ p refer to a complete set of quantum dot states, and W mfl is the coupling matrix element 
between channel m and the dot state [i. The inverse propagator D has the form 

D, U (E) = {E- EJS,,, W m ,W^, (4.32) 

m 

where the sum runs over m = L,R and where E^ includes an energy shift resulting from 
the coupling to the leads. First consider the system at zero magnetic field. Then the matrix 
elements can be chosen real. To simplify the argument, we restrict ourselves to two levels 1, 
2 in the dot. The quantum dot propagator is given by 

i / v 1 f E-E 2 + iT 22 /2 -iT 12 /2 , 

D (E) = — 1 1 , 4.33 

v ; detD(E) \ -iT 2l /2 E - E x + iF u 2 1 



with r ^ = 2tt J2m WmpW^u,. Combining Eqs. fl4.31|) and ( |4.33| ), one obtains the transmission 
amplitude t LR (E) between left and right lead 

2m 

tLR(E) = - deW{E) hn(E - E 2 ) + l22 (E - E{) + 

-(7n r 22 + 722rn - 7i 2 r 2 i - 72ir 12 )], (4.34) 

where 7^ = W^W^. Substituting the definit ions of and r y IMU into Eq. ( |4.34| ) , one finds 
that the term in the round brackets on the right hand side vanishes identically. Therefore, 
if 711 and 722 have the same sign, the transmission amplitude tin vanishes for some energy 
between E\ and E 2 . We thus obtain a transmission zero. Note, that for the above argument, 
it is crucial to keep the off-diagonal elements of the propagator. They were neglected in 
Ref. 131. 



It turns out that both the condition of time-reversal invariance and of single-channel 
leads are necessary for a transmission zero. If one of these conditions is relaxed, one generally 
no longer finds a vanishing transmission. E.g. consider the case that time-reversal invariance 
is broken. The matrix elements W mfl then can take complex values and the term in round 
brackets in Eq. ( f4.34|) is no longer guaranteed to vanish. As a consequence, one generically 
finds a smooth phase slip and not an abrupt phase jump between resonances. However, 
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the phase slip may take place on an energy scale much smaller that the total width T. 
This is suggested by the following observation. In the experiment, all resonances behave 
similarly (they display nearly the same phase and peak height). This suggests that different 
resonances have similar coupling matrix elements. For the simple two-level dot studied here, 
Jfiu ~ 7 and Y^ u « Y independent of the level indices /x, v. This choice results in a near 
cancelation of the terms in the round bracket in Eq. (|4.34p , and the phase slip may get very 
sharp. 

We finally discuss the effect of electron-electron interactions and finite temperature on 
the transmission zeros. Lee [IlJ argued that the analysis presented for the noninteracting 
case applies equally to interacting systems provided quasiparticle excitations remain well 
defined at E = Ep. The argument is based on the Friedel sum rule which is known to hold 
even in the presence of interactions. Unfortunately, no numerical or analytical study has 
explicitly demonstrated transmission zeros in the presence of interactions. Finite tempera- 
ture smears out the transmission zeros. The sharp phase jump is then replaced by a rapid 
but continuous lapse of the phase WJ\. The precise energy scale of this lapse depends on the 
detailed electron dynamics in the dot. We note that the phase lapse observed in Ref. is 
indeed continuous and has been resolved experimentally, cf. Sec. [Hj. 



3. Disordered dot 

The dots used in the experiments |||7] were roughly 50 times smaller than the elastic 
mean free path. Under these conditions, transport through the dots is ballistic. Most of the 
theoretical studies of the transmission phase pertain to the ballistic regime, and neglected 



the influence of disorder. Baltin and Gefen [32,48] took a different approach and investigated 



the transmission amplitude through a disordered quantum dot. They identified a generic 
mechanism for phase correlations and formulated an approximate sum rule: According to 
this rule, the change in the transmission phase AO between two consecutive transmission 
valleys is (mod 2tc). The frequency of deviations from this rule is small in A/Uq where 
A is the mean single particle level spacing and Uq the charging energy of the quantum dot. 
The sum rule pertains to individual, disorder specific systems. Baltin and Gefen also calcu- 
lated a disorder-averaged phase-phase correlation function and observed an enhancement 
of correlations with increasing interaction strength. Below we present the main arguments 
of the study @. 

The analysis of Baltin and Gefen proceeds from the scattering theory in the weak- 
coupling limit as presented in Sec. |1V A3|. However, the basic argument can be understood 



from a simple toy model |32| based on two assumptions: (i) The single-particle level spacing 
is constant S\ = AA and (ii) the product of the coupling matrix elements is a random variable 
Vj^Vj^* = Vrjx where f]\ takes the values +1 and —1 with equal probability. These simple 
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FIG. 16. Re t (left) and phase (right) for a specific sequence of resonances, for Uq 
panel) and Uq = 60 A (lower panel). Taken from Ref. |3^| . 
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assumptions model the fluctuations in the wave functions due to disorder. A more elaborate 
random-matrix calculation of the energies and the couplings confirms the predictions of 
the toy model described below. Consider first the noninteracting case (Uq = 0). The 
transmission amplitude then becomes 

Vx 



t(E) = vY: 



E-£ x + iT> 



(4.35) 



where the system can be tuned in or out of resonance by the gate voltage V g (we assume 
£\ = £f ] -eV g ). At each resonance, the transmission phase increases by ir. The signs of rjx 
govern the phase evolution between the resonances. If rj x ■ ?7a+i > there is a phase lapse by 
Ti between the resonances A and A + 1; for 77^ ■ ?7a+i < there is no phase lapse. In strictly 
one-dimensional systems rjx alternate in sign implying no phase lapse as discussed in Sec. 



IV B. For a disordered noninteracting dot the phase lapses occur at random. 



The situation is different in the presence of interactions. In the conductance valleys a 
large number of random terms ~ Uo/A contribute to t(E). The terms in the iVth valley 
and the (N + l)th valley differ very little from each other (essentially by the contribution 
of one level). The transmission in the valleys is therefore determined by a background that 
varies little between neighboring valleys. Figure [16] shows the evolution of Re(t) and arg(t) 
for for kT = A/12 and specific series of couplings, both for the noninteracting case Uq — 
(upper panel) and Uq = 60 A (lower panel). For interacting electrons the change in phase 
Aarg(t) = as expected from the sum rule. Note, however, that the phase evolution as a 
function of voltage is different from what has been seen experimentally: The phase in the 
disordered dot decreases near the resonances and not in the valleys between resonances. 
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FIG. 17. Correlation function Cg vs. valley number for kT = A/12. Interacting electrons 
(boxes and diamonds, Uq = 50A, x = x = 0.5) show strong correlations, non-interacting electrons 
(circles, x = x = 0.25) do not. Toy model statistics (boxes) agrees well with random matrix model 
statistics (diamonds). Inset: Increasing the interaction increases the correlations (upper curve: 
Uq = 50A, middle: Uq = 25 A, lower: Uq = 12 A). Taken from Ref. ||. 

The crucial influence of the interaction is further illustrated by the disorder-averaged 
phase correlation function 



C e = (cos6(x,N) cos 6{x, N + 8N)). 



(4.36) 



Here, x and x with < x, x < 1 specify the voltage in the N-th and (N + 5N)-ih valley. 
Cq decays slowly on a scale 5N ~ Uq/ A. Figure |l^ shows Cg vs. SN for kT = A/12. While 



there are no correlations for noninteracting electrons, one observes a slow decay of Cg in the 
presence of interactions. The decay is slower for stronger interactions as expected from the 
arguments given above. 

The sum rule found for disordered dots is in accord with the observations made in the 
experiment 0, however, the evolution of the phase vs. gate voltage is in disagreement. We 
believe that the latter fact indicates a fundamental difference between the phase behavior 
found in disordered (or fully chaotic) systems and the phase behavior in ballistic dots with 
mixed dynamics. 

D. In phase resonances 

Several mechanisms for peak-correlations and in-phase resonances have been proposed. 
These include the effect of finite temperature, deformations of the dot confining potential, 
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and energy shifts caused by the coupling between the dot and the leads. We summarize 
these mechanisms in Sees. fV D 1| - [IV D 4| . The mechanisms rely on the assumption that the 
transmission is dominated by a few quantum states of the dot. The dynamical origin for 
this behavior is discussed in Sees. [IV D 3[ - [TV D 5| . The investigation suggests that the in- 
phase resonances result from features of the short-time dynamics in the dot that completely 
dominate the quantum transport in the tunneling regime. Alternative proposals for peak 
correlations are reviewed in Sec. II V D 6. 



1. Finite temperature 

Transport at finite temperatures differs qualitatively from transport at zero temperature. 
Close to a conductance peak electrons tunnel through a single level if kT <C A. In contrast, 
few or many levels contribute to each conductance peak if kT is of the order A. Similar sets 
of levels contribute to neighboring peaks which generates peak correlations. Oreg and Gefen 
f49|| studied the transmission phase correlations resulting from this effect. More recently, 



temperature induced correlations of the conductance peak-height in fully chaotic or diffusive 



dots have been calculated [5(| using random matrix theory. A comparison with the peak- 



height correlations measured for chaotic dots has been given in Ref. BT 



Oreg and Gefen demonstrated in-phase behavior for a simple model with two levels in 
the dot. The calculation is an application of the formalism developed in Sec. [IV A| . Due to 



the small number of levels, all contributions to the transmission amplitude can be computed 
explicitly. The starting point is the two-level Hamiltonian for the quantum dot, 

#qd = £ a a j a + £ b b ] b + U a ] ab ] b, (4.37) 

where Uq denotes the strength of the electron-electron interaction. The dot is coupled via 
matrix elements V^, V a ~, V b L , V b R to two single-channel leads L and R. As in Sec. |IV A 



all energies are counted from the Fermi energy in the leads. To obtain the transmission 
amplitude (|4.20|) through the dot, we compute the occupation probabilities of the dot states. 



The quantum dot supports four many-body states. They can be labeled |0) (no electron in 
the dot), | a) = a^|0), \b) = b^O), and \ab) = a^b^O). The respective energies are E = 0, 
E a = £ a , Eb = £b, and E a b = £ a + £b + U$. The equilibrium probabilities to find the dot in 
either of the states \i) = |0), \a), \b), \ab) are given by 

P, t = exp(-PE i )/ J Eexp(-PE i ). (4.38) 

i 

In terms of the Pi, the matrix elements of the retarded Green function are 

G aa (E) = (Po + Pa) = * „ /9 
E - £ a + iT a /2 
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FIG. 18. Left: Occupation probability of the four many-body states as a function of V g . Right: 
Argument of tRL (solid line) and its magnitude (dashed line) as a function of V g . For the parameters 
A = lOfieV, kT = 20fieV, U = 500/xeF. Taken from Ref. Q. 
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(4.40) 



where T a:b = \V a L b \ 2 + \ V a R b \ 2 . Substitution of Eqs. ( [Op , ( piOl) into Eq. yields the 

transmission amplitude £qd- 

Two consequtive resonances with similar phase can arise in this model if one level is 
significantly stronger coupled to the leads than the other level. For definiteness, let a be the 
strongly coupled level. At the first resonance, either an electron or a hole tunnel via level 
a through the dot while level b is empty. These processes are represented by the first term 
on the right hand side of Eq. ( 4.39|) . At the second resonance, again charge is transfered 
predominantly via level a, this time with level b being occupied. The corresponding processes 
are described by the second term on the right hand side of Eq. (|4.39|) . Two consecutive 



resonances are then dominated by the same single-particle level a. This scenario differs 
qualitatively from the transmission of independent electrons (which is recovered in the limit 
Uq — > 0), where two consecutive resonances are associated with two different single-particle 
states. 

We use the parameterization £ a = V g , £ b = V g + A of the single-particle energies in 
terms of the voltage V g at the dot. The occupation probabilities of the four many body 



states and the modulus and phase of £qd are displayed vs. V g in Fig. [18| For both plots, 
A = 0.5kT, and \V b L \ = \V b R \ = |0.3V^ L | = |0.3V^|. The two resonances are in phase since 
both are dominated by the same level a. The phase lapse between resonances takes place 
on the scale T a . This is expected since we neglected the off-diagonal elements of the dot 
Green function, cf. Sec. [IV C| . 
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FIG. 19. Energy levels of the two-dimensional harmonic oscillator with frequencies u x = to 
and uj y = u>(l — (3) as functions of 0. Taken from Ref. [52|. 



2. Deformations 

Transport experiments in the Coulomb blockade regime study the transmission through 
quantum dots as a function of a gate voltage V g applied to the dot. Most theoretical studies 
of the Coulomb blockade assumed that V g exclusively regulates the depth of potential well in 
the dot. Changes in V g then translate into an overall shift of the single-particle levels of the 
dot. This picture is clearly an oversimplification. In any experimental realization of a quan- 
tum dot, variation of V g modifies both the depth and the shape of the dot confining potential. 
The consequences for the Coulomb blockade were first investigated by Hackenbroich, Heiss, 
and Weidenmuller Shape variations lead to avoided crossings of single-particle levels. 



This can generate sequences of conductance resonances which carry essentially the same 
internal wave function. Within such sequences one may find in-phase resonances and/or 
strong correlations of the conductance peak heights. Vallejos, Lewenkopf, and Mucciolo ]53]| 
showed that shape deformations may change the spacing distribution of Coulomb blockade 
peaks. Stopa []54j] demonstrated shape deformations in self-consistent calculations of the dot 
confining potential. We discuss shape deformations for a quantum dot with a parabolic con- 
fining potential. Then we consider more realistic confinement potentials with level repulsion 
and chaotic classical electron motion. 

Consider a quantum dot with a parabolic confining potential. The potential is charac- 
terized by two oscillator frequencies uj x and u y . For simplicity, we keep u x = ui fixed while 



uj y = co(l — 0) depends on the deformation parameter (3 with < /3 < 1. Fig. [19] shows 
part of the single-particle spectrum versus (3. The levels form a network of intersecting 
straight lines. The single-particle state associate with each line is characterized by two 
non-negative integer quantum numbers (X x , X y ). States with small (large) values of X y have 
small (large) negative slopes. We show below that a pattern similar to Fig. ITU] may also 
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Pn Pn+i Pn+2 

FIG. 20. Schematic illustration of the mechanism described in the text. At the deformations 

Pn, Pn+i, Pn+2 one additional level (indicated by the arrows) is occupied. Taken from Ref. |j52[ , 

be found for the spectrum of deformed non-integrable potentials, with one difference: The 
points of intersection disappear and are replaced by avoided crossings. 

To illustrate the mechanism for peak correlations, we replace Fig. [19] by the idealized 
picture shown in Fig. A set of equally spaced straight lines (A-levels) runs nearly parallel 
to the /5-axis, a second set (5-levels) has large negative slope. To model the generic case, 
actual crossings have been replaced by avoided crossings. The distance between _B-lines is 
denoted by 5j3 CTOSS . Wave functions retain their identity across avoided crossings. Therefore, 
the wave function on any A-level is nearly independent of (3. A change of the gate potential 
and, therefore, of (3 is slow on the scale of the characteristic times of the quantum dot. 
Therefore, each electron in an occupied level follows the deformation (3 adiabatically. At 
each avoided crossing, the electron wave function changes from A-type to B-type or vice 
versa. 

Suppose that for some (3 the N lowest single-particle levels of the dot are occupied and 
that the highest occupied level is an A-level. What happens as (3 is increased by <5/3 croS s? 
The last occupied level moves adiabatically down at the avoided crossing. The associated 
wave function switches from A-type to B-type. The A-level that was the last occupied 
level before the avoided crossing becomes empty. Upon further increase of f3, this very same 
A-level can be occupied by another electron (at deformation Pn+i)- This process can repeat 
itself for a second, third, etc. time ifSP cross ~ 6P e , the change of P needed to pull another 
electron into the dot. Then, subsequent conductance peaks would not be independent, but 
rather be manifestations of essentially the same single-particle state. Strong correlations of 
the conductance peak heights and the transmission phases are expected for such sequences. 



How does the idealized picture of Figs. 20 change for a more realistic spectrum with 



level repulsion? As a generic example consider the Hamiltonian 
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FIG. 21. 
Ref. p|. 



P 

Energy levels as a function of of the Hamiltonian ( [4.41| ) for p 



= 0.004. Taken from 



H = ^ + ^f(x 2 + (1 - /3)V) - P huL\ (4.41) 

where L is the dimensionless ^-component of the angular momentum operator. The three- 
dimensional analogue of H is known as the Nilsson model and has been quite successful in 



explaining the spectra of deformed nuclei |55| . For f3 > and p ^ 0, H is not integrable 
and displays level repulsion. 

Figure |2l] shows part of the spectrum of H versus (3. The overall pattern is quite similar 
to Fig. [19|, however all crossings have been replaced by avoided crossings (mostly not visible 
in Fig. |2l"l) . One observes a set of nearly flat lines. Numerical inspection shows that the 
associated wave functions retain their identity over a large range of (3. The intensity of these 
wave functions is concentrated along periodic orbits which oscillate in the x-direction, with 
little or no motion in the ^/-direction. One expects that the bulk part of the current through 
the dot is carried by such orbits. The numerical observation of stable periodic orbits in the 
Nilsson model of Ref. [52] provided the first evidence for the dynamical origin of strongly 
coupled quantum states. We note that Fig. ^l] also reveals steep levels. The corresponding 
wave functions change strongly with deformation as expected for parametric variations in 
chaotic systems. 



3. Integrable dot 

In the preceding sections we presented qualitative arguments for peak correlations in- 
duced either by temperature or deformations. A synthesis of both approaches was given by 
Baltin et al. pHJ . It was shown that a quantum dot of the shape of a deformed harmonic os- 
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cillator can support sequences of up to 30 conductance resonances with similar transmission 
phase and similar peak height. All resonances within such a sequence are dominated by a 
single strongly coupled eigenstate of the dot. The study ]3TJ] pertains to an integrable dot 
in the regime of ballistic transport and is therefore restricted in generality and universal- 
ity. However, aspects of the study are also relevant for ballistic dots of more general shape 
provided these dots support short orbits that strongly couple to the leads. 

Consider the tunneling Hamiltonian H defined in Sec. IV A and assume that the dot 
confining potential is an anisotropic harmonic oscillator potential. The states of the dot can 
be labeled by two quantum numbers A^, X y and the energy eigenvalues E\ for A = (X x , X y ) 
are given by 



£ x = = hu x {X x + i) + ruv y {V g ){X y + i) - aV g + Eq. (4.42) 



To describe the deformation, we assume that the oscillator frequency u) y (V g ) = U) x (l — j(V g — 
Vq)) in the transverse direction y depends linearly on the gate voltage V g while the frequency 
uo x in the longitudinal x-direction is held fixed. The parameter a relates the overall depth 
of the dot potential to the gate voltage. The constants Eq and Vq determine the number of 
electrons on the dot at zero deformation. 

The matrix elements V k L x , V kX for tunneling from the left and right lead to the quantum 
dot are given [56] by the integrals 



2 



d<f>x x ,xJx,y) dip k (x,yY 



dx x ' v dx 



(4.43) 

J X=X B 



I ( Ft) 

where ip k denotes the wave function with wave vector k in the left (right) lead , and 
<&\ x ,\ the wave function in the dot. The integration extends in the y-direction and xb is 
arbitrary but must be located within the barrier. We restrict ourselves to the case of a single 
transverse channel in each lead. The nodes of the wave functions of flat (steep) levels with 
large X x (X y ) are predominantly carried by the x-component (y-component, respectively). 
Thus, the wave functions of flat levels extend much further into the barrier region and have 
considerably larger matrix elements V k ^ than those of the steep levels. This important 



property is illustrated in Fig. |22|. Due to their stronger coupling, flat levels will carry larger 
current that steep levels. 

According to the scenario discussed in Sec. |1V D 2| a flat level of the dot gives rise to 
a sequence of Coulomb peaks when it undergoes a single level crossing between every two 
successive peaks. At finite temperature, this condition need not to be met exactly, but must 
hold on average for a sufficiently large number of peaks. Within our model, we can estimate 
the number AN of correlated resonances. For Uq ^> A, the distance between neighboring 
resonances is SV g = U /a. The number of intersection points of a flat level (A x . , X y = 0) 
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FIG. 22. The thin solid line shows a cross section of the potential in longitudinal direction, the 
two barriers lying at opposite ends. The overlap of the dot wave functions (probability shown as 
dashed lines) with the lead wave function (probability shown as a solid line on the left) increases 



strongly with the quantum number X x . Taken from Ref. [31 



with steep levels is found using Eq. ( [4.421) . For the total number of crossings in the voltage 
range V g — V one finds 



K l{V 9 -Vo) 
2 1 - 7 (V S - V 



(4.44) 



The number of steep levels around a given energy increases with with deformation. This 
yields an increase of N c , and causes the divergence for the unphysical case of extreme de- 
formation ujy ~ 1 — j(V g — Vo) — > 0. One crossing within the interval SV g occurs for 
(dN c )/ (dV g ) = a/Uo- This condition yields the voltage V* where maximal correlations of 
the Coulomb peaks are found; V* is used below. The number AiV of correlated peaks can 
be estimated as the number of resonances for which the flat level (X x , 0) stays within the 
energy interval A around Ep. This yields 



AiV 



a 



2jU 



(4.45) 



with the deformation u y /u x = \J{jUo)/{2a) X x . The result ( [4.45Q explicitly relates the 
number of correlated levels to our model parameters. We note that AiV sets an upper bound 
for the number of correlated peaks. If V g is substantially different from V*, the number of 
level crossings no longer matches the number of Coulomb peaks. Then the sequences of 
correlated peaks are shorter than AiV and the correlations are weaker. 

The two-terminal conductance G and the transmission phase through the quantum dot 



are calculated using the formulas derived in Sec. [IV A 3| . In Fig. |23] we show G as a function 
of the gate voltage V g for two values of kT. The results are obtained for a = 1, 7 = 0.005, 
E = —11 and Vq = 90 where the energies and voltages are measured in units of Uq. The 
single-particle level spacing is A = 0.03£/ which roughly corresponds to the experiments 
30|,[7[|. About 100 Coulomb blockade resonances occur in the interval 100 < V g < 200. In 
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FIG. 23. Conductance G vs. gate voltage for (a) kT = A/5 (left) and (b) kT = 2 A/5 (right). 



Taken from Ref. [31]. 



both plots strong peaks with similar peak heights appear whenever a flat level is close to 
the Fermi energy. The peak-height correlations are more pronounced at higher temperature 
(case (b)). In the regions V g < 130 and V g > 160, there is no the flat level close to the 



Fermi energy and G is small. On the scale of Fig. 23, some of the conductance peaks are 
not visible. 

The phase 9 of the transmission amplitude is shown vs. V g in Fig. |24]. For the parameters 
chosen, the transmission is dominated by the flat level A = (14,0); the width of this level 
is = A/15. The solid lines at the bottom of the plot show the conductance peaks and 
help to identify the resonance positions. One observes a strikingly similar behavior of the 
phase at all resonances. This behavior is found not only within the V g interval shown but for 
the entire interval 130 < V g < 160 comprising 30 resonances. The phase increases by 7r at 
each resonance and displays a sharp lapse by 7r between adjacent resonances. The increase 
at resonance occurs on the scale kT (we assumed kT > T) and the phase lapse between 
resonances on the scale T (see the discussion in Sec. |IV C|) . 

The identical behavior of 9 at all resonances reflects the fact that at each resonance, 
the transmission through the dot is dominated by a strongly coupled level A. Similarly, 
the phase lapse between adjacent resonances is caused by the dominant level: At finite 
temperature the level A has a finite probability of being either occupied or empty and, thus, 
contributes to both an electron-like and a hole-like cotunneling process. The contribution 
of both processes to the transmission amplitude is 
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FIG. 24. Phase 9 of the transmission amplitude versus gate voltage V g at kT = A/5. The solid 
lines at the bottom of the plots display the conductance peaks. The flat level X x = 14, Ay = is at 



or near the Fermi energy. Taken from Ref. [31 
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E-(£ x + (N-l)-U )+iT x /2 



(4.46) 



where the first (second) term represents the electron (hole) contribution, respectively. As 
the gate voltage V g scans the N th conductance valley, the sign of Re(tQD) reverses, leading 
to a lapse in the transmission phase. 

We emphasize that the phase lapse between resonances is a genuine interaction effect. 
For vanishing charging energy, the cotunneling amplitude ( f4.46|) would reduce to a single, 
temperature-independent term. The phases of the transmission amplitude in consecutive 
valleys would not be correlated, and there would be no systematic phase lapse between 
resonances. We also note that the systematic phase slip only occurs at finite temperature. 
At zero temperature a flat level contributes to either particle-like or to hole-like cotunneling 
and no phase lapse is expected within the present model. 



4- Energy shift 



Silvestrov and Imry |58[ demonstrated a mechanism for in-phase behavior that requires 
neither temperature nor dot deformations. They showed that the energy shift resulting from 
the coupling between the dot and the reservoirs plays an important role and may cause in- 
phase resonances. Energy shifts become relevant when they are of the order of the mean level 
spacing. This requires a coupling of the order T \n(U /T) ~ 2nA. As a result, the scenario of 
Silvestrov and Imry cannot explain the phase correlation in the Yacoby-experiment (where 
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FIG. 25. Projection of the scattering wave function corresponding to a broad level for three 
different energies. The energies are indicated by the numbers in the upper left corners. The broad 
level is the superposition of two contributions which are quantized on different classical trajectories. 
One trajectory (left plot) has its classical turning point at the left contact, yielding strong coupling 
to the lead. Taken from Ref. [pal . 



r <C A) but is possibly relevant for the strong-coupling regime explored in the Schuster- 
experiment. The number of peaks that may be correlated due to energy shifts is estimated 
to be ~ r/(2vrA) \n(U /T). 

Slivestrov and Imry illustrated their idea for a quantum dot of linear dimension I and 
(dimensionless) confining potential 



V(x,y) = -Ax 1 




(4.47) 



The dynamics in the dot is non-integrable but the potential is approximately symmetric 
similarly to the dot used in the experiments. The dot is connected to a single lead attached at 
x = and extending in negative x-direction with a parabolic confining potential in transverse 
direction. The Schrodinger equation in the dot and the lead can be solved numerically on a 
lattice. The kinetic energy is modeled by a standard nearest neighbor hopping term. Within 
the energy interval 1.5 < e < 4.7 used below only a single propagating mode exists in the 
lead. The energy range corresponds to above-barrier scattering through the dot. 

The numerical solution of the single-particle scattering problem reveals scattering res- 
onances with very different width. Silvestrov and Imry extracted the widths by fitting the 
numerical results with a sum of Breit-Wigner resonances. This yields broad resonances with 
T of order (5 — 7) A as well as very narrow resonances with width much less that A. The 
projection of the scattering wave function corresponding to one broad resonance is shown 
in Fig. It has large probability near the contact to the lead, and may be understood as 
the superposition of two contributions each quantized on a short classical trajectory. 

Motivated by their numerical findings for the single-particle problem, Silvestrov and 
Imry then investigated many-particle effects. They studied a model with only one level 
(denoted //) coupled strongly to the lead. The width of this level was assumed to be T ^> A, 
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the width of all other levels is much smaller than A. Even though one level was strongly 
transmitting, the standard charging model was employed to model the electron-electron 
interactions. It was assumed that the charging energy Uq ^> T. To understand the origin 
of the energy shift calculated below, it is useful to decompose the total occupation number 
N = N' + in two part, the first part N' = ^)w„ fi\ describing the total occupation of the 
sharp levels and the part the occupation of the broad level. The charging energy then 
reads 

U = ^U (N 2 — N) — ^U (N' 2 - N') + UoN'h^. (4.48) 

For fixed N', the charging energy contributes a term ~ n M to the Hamiltonian. This term can 
be combined with the single-particle term S^h^ and yields the effective energy S^ N '^ = £ M + 
UqN' which depends on the occupation number N'. Here and below we use the superscript to 
indicate the number of electrons in the sharp levels. For fixed N' and vanishing coupling of 
the sharp levels, the scattering problem through the dot reduces to the well-known problem 
of a single state interacting with a continuum. The exact solution of this problem is known 



|55| , |58[| . The coupling between the state and the continuum reduces the total energy of the 
system below its value found for the decoupled case. The energy shift introduced by the 
coupling is also known exactly |55l j58|1 . 

We now investigate the ground state of the dot as a function of the gate voltage. Assume 
that the dot is initially filled with N electrons and that the gate voltage is then tuned beyond 
the charge degeneracy point for N and N+1 electrons. The states in the dot are labeled in 
the order of increasing energy, and it is assumed that // > N + 1. Neglecting the coupling 
to the continuum, the ground state of the dot with N + 1 electrons is obtained by filling the 
lowest N + 1 levels and leaving the broad level \x empty. When the coupling is included, the 
total energy of this configuration is given by 

N + 1 1 p / AT? \ 

= E £a + - 2 U [(N + if - (N + 1)] - A- In J , (4.49) 

where we have omitted the energy of the electrons on the leads (energies are counted from 
the Fermi level). The last term on the right hand side is the energy shift resulting from the 
coupling between the dot and the lead. The expression is valid provided £^ N+1 ^ 3> T. The 
corrections for |^ Ar+1 - ) | < T can be found in Ref. [5S|. The energy ( |4.49| ) may be compared 
with the energy 

v 1 ^ UEf 



4? = E £a ^U [(N + l) 2 - (N + 1)] - - In [ 



(TV) 



(4.50) 



of the configuration where the broad level and N' = N sharp levels are occupied. In the 
weak-coupling limit, E^ +1 ' is always less than E^ t \ and the transmission at the N+1 
Coulomb peak proceed via level N + 1. However, when T is large enough so that 
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" £n+i < In , (4.51) 



the configuration for which level /x is occupied has lowest energy. In this case, one finds a 
broad Coulomb peak caused by the transmission through the level fi. 

When the voltage is further increased, the two functions E^P and E^ t +1 ^ cross. The 
voltage at the crossing point is determined by the equation 

£,{V g ) = U N Uo tt^—. (4.52) 

exp[2iT(£^ - £ N+1 )/r\ + 1 

At the crossing point, the ground state jumps onto the branch E^ and the current- 
transmitting level \i is again empty. At zero temperature, the jump from one electron 
configuration to the other is accompanied by a sharp jump by tt of the transmission phase. 
The process of filling the broad level and jumping to another electron configuration can 
repeat itself. The number of consecutive correlated resonances caused by this process is 
~ r/(27rA) ln(Z7o/r). We note that the analysis described above pertains to the ground state 
of the quantum dot. A generalization to finite temperature is required for a quantitative 
comparison with the experiment (e.g. to explain the finite scale observed for the phase lapse). 



5. Bouncing-ball tunneling 

In Sees. [IV D 2I4TV D4| we showed that quantum dots with integrable or mixed classical 
dynamics may support a subset of quantum states with exceptionally strong coupling to 
the leads. These states are quantized on short periodic orbits that connect the contacts to 
the leads. The dynamical origin for strong coupling was further clarified by Hackenbroich 
and Mendez ||59|| . Using a tunneling Hamiltonian approach, they showed that the relative 
coupling strength of the dot states strongly depends on the transverse width of the tunneling 
region between the dot and the leads. Exceptionally strong coupling to few dot states is only 
found for sufficiently wide leads, ka e s > 1. Here, k is the Fermi wave number and a e g the 
effective transverse width of the tunneling barrier. The transport under these conditions may 
be termed bouncing-ball tunneling (BBT), as the strongly conducting states are quantized 
on classical trajectories bouncing between the contacts to the leads. A unique fingerprint 
of BBT is found in the regime of strong coupling V > kT: Then the peaks in the tail 
of a sequence of correlated Coulomb peaks develop a characteristic line-shape asymmetry. 
The origin of this asymmetry is the breaking of the particle-hole symmetry as the strongly 
coupled bouncing-ball state moves away from the Fermi energy. 

The crucial role of the transverse barrier width a e g can be understood from the following 
simple argument: Confinement in a lead of width a e g yields the transverse momentum 
spread h/a e s for electrons injected in the quantum dot. Wide leads therefore result in near 
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normal injection and provide exceptionally strong coupling to the bouncing-ball states. In 
the tunneling Hamiltonian, a e g enters in the the matrix element for tunneling between 
the lead / and the state A in the quantum dot. For sufficiently high barriers, Vjjj is given by 



h 2 



dsi/ji(s,z)d z ip* x (s,z) 



2=0 



(4.53) 



where the integration is performed along the edge between the potential barrier and the 
quantum dot (d z denotes the derivative normal to the barrier). The wave function ip\ 
corresponds to Dirichlet boundary conditions in the dot, while the barrier tunneling is fully 
included in the lead wave function ipi. The transverse potential in the tunneling region 



can be taken quadratic |JD| yielding ipi ~ cjexp[— (s — si) 2 /2a 2 s ], where s is the transverse 



coordinate, si the center of the constriction and a e g its effective width. One can restrict the 
calculation to the lowest transverse mode since higher modes are suppressed by the barrier 
penetration factor (included in q). 
The central result of Ref. 



<fc = r?r?/(r 



"a 



rf: 



j9| is presented in Fig. |26|. Shown is the coupling strength 
over a sequence of 100 quantum states labeled by the index A. Note 



that g\ is proportional to the conductance peak height G\ = (e 2 /h){-K /2kT)g\ measured 
in low temperature Coulomb blockade experiments The dot is described by a hard- 
wall confining potential at the boundary, parameterized in polar coordinates by R(<p) = 
Rq[1 + ecos(20)]. Here, e measures the quadrupolar deformation out of circular shape. A 
nonzero value for e mimics the dots used in Refs. ||,[7j]. As in the experiments, the leads 
are attached opposite to each other at the boundary points closest to the origin (the points 
with = ±7r/2). For the value e = 0.2 used for Fig. ^ the classical dynamics in the dot is 
almost completely chaotic except for two large islands associated with stable bouncing-ball 
motion between the contacts to the leads. 



For narrow leads (Fig. p6(a)) one observes many peaks with widely different peak height. 
In striking contrast, the results for wide leads (Fig. ^(b)) show a few isolated large peaks, 
separated by 15 — 25 levels with much smaller peak height (not visible on the scale of Fig. 
|26](b)). All large peaks are associated with states quantized on stable bouncing-ball orbits. 
This is illustrated in the inset for the state A = 337. The height of the small peaks not 
resolved in Fig. |2B|(b) is typically two or more orders of magnitude smaller than the maxi- 
mum peak height. Such tiny peaks are difficult to resolve in Coulomb blockade interference 
experiments. Interference experiments with wide leads are therefore only sensitive to the 
strongly coupled bouncing-ball modes. 

We now turn to the calculation of the transmission coefficient Tqd and the phase 6qb- 
We assume wide leads and kT < A. The case of a weakly coupled dot was studied in Sec. 
|1V A| . Here, we consider a more open dot characterized by ~ A. We identify fi with the 
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FIG. 26. Coupling strength g\ for a sequence of 100 states in the interval 290 < A < 390. 
Results are for a quantum dot with quadrupolar shape and deformation e = 0.2. (a) Narrow leads: 
ka eS = 0.1. (b) Wide leads: ka eS = 5.0. Here k = (2m*E x ) 1/2 /h is the wave number associated 
with the bouncing-ball state A = 337. Inset: Real-space projection of the state A = 337. 



bouncing-ball state closest to the Fermi energy. All other states A ^ fi in the vicinity of Ep 
have a much smaller width T\ ^ A. The Green function for this case may be obtained using 
the equations of motion method (see Sec. |1V A 2| ). It is diagonal up to small off-diagonal 
corrections 0(JT\T tl /A), and given by 



N'=0 



N> 



E-{8„ + U Q -N>)+iTj2 



(4.54) 



Here N' counts the total number of electrons in all dot levels except for the level /1 and 
Pn> is the respective occupation probability. The transmission through the states A 7^ ji 



is negligible. Equation ( 4.54|) is the generalization of the weak-coupling result derived in 



Sec. [IV A 2| to the case of a single strongly conducting quantum state. The probability 
P/v = Z~ l exp[— Q(N') / kT] with Z = J2n' exp[— Q(N')/kT] is related to the thermodynamic 
potential ft(N') of the dot. To evaluate P/y/ we replace fijv' by fl% + [£ M + Uq ■ N}(n^)i 
where QP N , is calculated for the dot with level \x excluded from the spectrum, and 



>N' 



N' 



dE Im 



f(E) 



E-^ + Uo-N^+iT^ 



(4.55) 



is the canonical occupation probability of level \x. 

In Fig. |27] we show the transmission coefficient Tqd and the transmission phase #qd vs. 
gate voltage V g (we assumed £ M = — eV g ). All peaks shown result from transmission 
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FIG. 27. (a) Transmission coefficient Tqd and (b) phase #qd as a function of gate voltage V g 
evaluated for kT = 0.2A, = 1.5A, and Uq = 12A. The peak asymmetry in (a) is a unique 
fingerprint of BBT. 



through the level /i = 337. The peaks have comparable height and similar phase in qualita- 
tive agreement with the experiments ||[7|. The central peak has a Lorentzian shape of width 
T^. Note that the transmission peaks develop a peculiar asymmetry as the conducting level 
moves away from the Fermi energy: Each peak to the left and to the right of the central 
peak falls off more rapidly on the side facing the central peak. This pattern extends over the 
whole sequence and becomes more pronounced for the peaks in the tails. The asymmetry 
results from the breaking of the particle-hole symmetry in the transmission through the 
bouncing-ball state and is a unique fingerprint of BBT. Inspection of the data of the exper- 
iment |7| reveals the same asymmetry, providing strong evidence that the peak correlations 
in this experiment are due to BBT. 

Note that the magnitude of the conductance peaks in Fig. p7)(a) decreases in the tails 
of the sequence. This is in contrast to the experiment [[?] where similar peak-heights were 
observed. This discrepancy may have several reasons: First, the billiard model used here 
completely neglects deformations of the dot shape. As discussed in Sec. |I V D 2| , |1VD 3 



deformations can enhance peak-height correlations by "pinning" the conducting levels close 
to the Fermi energy. The model further neglects the electrostatic influence of the plunger 
gate on the point contacts defining the tunneling barriers. In most experiment, this influence 
is significant: Increasing the plunger voltage, inevitably opens the dot. To model this effect 
one may reduce the barrier height with increasing voltage. In turn, this would result in 
enhanced peak heights on one side of a sequence. This voltage induced enhancement and 
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the suppression by temperature may result in a number of peaks with similar peak height. 



6. Other mechanisms 



In two recent papers Wu et al. and Kang proposed that the in-phase resonances 

reflect an interference effect in the AB ring rather than a property of the quantum dot. The 
authors investigate the transmission of non-interacting electrons through a single-channel 
AB ring. A two-terminal AB ring is studied in Ref. and both a two-terminal and a 
four-terminal ring in Ref. ||62|| . The quantum dot is modeled in both studies as a symmetric 
ID double barrier well (identical barriers on either side). Full coherence is assumed through 
the dot and the AB ring. For an integer value of flux quanta threading the ring, a set of 
conductance peaks is found upon variation of the dot potential. All peaks are in phase. A 
new set of in-phase peaks appears at half-integer flux, while the previously found peaks 
disappear. Each peak of the new set is located between two peaks of the previous set. The 
explanation of the effect is straightforward and follows from the scattering theory presented 
in Sec. [IIIB| : There it was shown that the oscillatory part Gab of the conductance close to 
a resonance of the quantum dot is of the form 



Gab oc 1 + cos(0 - x), 



(4.56) 



where <fr — 2n^/^ is the dimensional flux and where % takes the values and 7r for even and 
odd states, respectively, of the quantum dot (cf. Eq. ( |3.29|) ). As a result, only even states 
contribute to the AB current if an integer number of flux quanta are threading the ring. 
Clearly, all even states are in phase. Conversely, for half integer number of flux quanta, only 
odd states contribute and, again, all with the same phase. Based on Eq. ft4.56|) one expects 
contributions of both even and odd states for intermediate values of the flux. Precisely this 
has been reported by Kang |62| for $ = $o/4 and $ = 3$o/4- 



The ideas presented in Refs. fl6"l"l]6"2|1 cannot account for the experimentally observed 



in-phase behavior. First, Refs. |n],|62| find in-phase resonances only for integer or half- 
integer number of flux quanta threading the ring. For generic values of the flux, neighboring 
resonances are always out of phase by tt. Second, the destructive interference responsible 



for the in-phase resonances in |61|j62]] is only complete if the transport through the AB 
ring is fully coherent. In the experiments |||7[], however, most of the current is incoherent. 
The coherent component comprises only about 10 or 20 percent of the total current. The 
incoherent ring current displays resonant behavior regardless of the value of the external 
flux or the parity of the dot state. Hence, experimentally both even and odd resonances are 
observed for all values of flux, and each current resonance is associated with the charging 
of the dot by one additional electron. The prediction [0,0 of missing resonances and 
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the addition of the charge 2e between neighboring resonances is in contradiction to the 
experiments. 
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V. CONTROLLED DEPHASING 



In 1996 Gurvitz |?kJ realized that mesoscopic electron devices could be utilized to study 
and controll decoherence. He investigated the resonant tunneling through two capacitevly 
coupled quantum dots. Gurvitz showed that the second dot could serve as a detector for the 
charge accumulated in the first dot. Solving the equations of motion of the entire system 
and then tracing over the detector variables, Gurvitz obtained the reduced density matrix 
of the measured dot. The coupling to the detector dot lead to a damping of the off-diagonal 
elements and thus caused decoherence. 

Gurvitz' paper has triggered new experimental and theoretical developments exploring 
the controlled decoherence of quantum transport by the environment. Since the destruction 
of coherence is not necessarily related to energy relaxation, the new subject is commonly 
referred to as controlled dephasing. The first demonstrations of controlled dephasing were 
achieved in beautiful which-path experiments of Buks et al. || and Sprinzak et al. |§. The 
experiments employed the transport through a single || or through two |§ tunnel-coupled 
quantum dots. In close proximity and capacitevly coupled to one dot was a quantum point 
contact (QPC). Electrons passing through the quantum dot interact with electrons in the 
QPC. This modifies the transmission through the QPC, so that the QPC serves as a detector 
for the charge of the dot. It was found that a current flowing through the QPC leads to 
dephasing in the quantum dot. In the experiment || dephasing was detected as a reduction 
of AB oscillations using a device with the quantum dot embedded in an AB ring. The set- 
up can be viewed as a which-path interferometer, the archetype of a position measurement 
apparatus in a double-slit device. We discuss the which-path experiments in Sec. [V A| . The 
theoretical work devoted to dephasing in mesoscopic electron structures is summarized in 
Sec. IVT3L 



A. Which path experiments 

1. Experiment of Buks et al. 



Fig. [28| shows the scanning electron micrograph and a schematic description of the set- 
up used in the which-path experiment of Buks et al. || . Part of the device, including the 
multi-terminal AB ring and the quantum dot embedded in its right arm is identical with 
the setup used in the Schuster-experiment (see Sec. |I]). The new element is the QPC to 
the right of the quantum dot. It is visible as a small constriction in the two-dimensional 
electron gas. The width of the QPC and hence its transmission coefficient is controlled by 
the gate voltage V g . The central gate is contacted with a metallic air bridge. This gate 
at the potential Vp depletes the area between the two arms of the ring and, at the same 



58 



time, serves as plunger gate of the quantum dot. The dot is tuned to the Coulomb blockade 
regime with the resistance of the tunneling barriers at either side of the dot being greater 
than h/2e 2 . 

The collector current J c measured for a fixed a.c. emitter voltage V# = 10// V displays AB 
oscillations as depicted in Fig. 29|. The solid line shows the AB oscillations with vanishing 
drain-source voltage Vd = across the QPC. The dotted line is measured with Vd = 100//V. 
The reduction of the visibility reflects the dephasing introduced by the which-path mea- 
surement . 

To study the dephasing quantitatively as a function of the QPC characteristics, Buks et 
al. investigated the calibration device shown in the inset of Fig. ^(a). The device contains 
a QD and a QPC similar to these in the which-path interferometer. The conductance of the 
QD is scanned through a series of Coulomb peaks by changing the plunger voltage Vp. Due to 
the proximity of QD and QPC, the transmission coefficient Td of the QPC is also affected by 
the potential of the QD. This gives rise to the smooth increase of Td with increasing plunger 
voltage Vp (see Fig. |50|(a)). However, whenever a conductance peak is being scanned and an 
additional electron is being added to the QD, Td displays a faster and opposite change with 
amplitude ATd on the scale of the peak width. This decrease of transmission reflects the 
change in the QD potential due to the additional electron. Fig. |30](b) shows AT d averaged 
over several Coulomb peaks as a function of the transmission 7^. The reduced value of ATd 
near 7^ = and % = 1 is a consequence of approaching the conductance plateaus. At the 
plateaus, changes in the QPC potential have very little effect on the QPC transmission. 

Fig. |31|(b) summarizes the main result of the experiment. The visibility of the AB 
oscillations is shown vs. the QPC voltage V g . The visibility is obtained upon dividing the 
peak-to-valley value of Ic by the average Ic- The measurement was done when the QD was 
tuned to a conductance peak using the central metal island as a plunger gate. Fig. j3l|(a) 
shows that the detector transmission increases from to 1 as V g is being changed. We note 
that any modulation of the visibility reflects the dephasing introduced by the detector; a 
small visibility corresponds to strong dephasing. For small detector voltage Vd = 10/jV, 
the visibility is practically constant with V g . For Vd = 100// V, one observes a pronounced 
structure: The visibility peaks near the conductance plateaus Td = 0, 1, and in between for 
Td = 0.5. It turns out that the weak dephasing at these points is due the small amount 



of which-path information in the detector, cf. Sec. [VB 1| . The detector is inefficient both 
at the conductance plateaus and at Td = 0.5 where the QPC current is very noisy. The 



theoretical approaches summarized in Sec. [VB| obtain for the visibility the result v = UQUd, 
where vq is the intrinsic visibility of the AB ring and v& the reduction of the visibility due 
to the detector. At zero temperature, v d is calculated to 

eV d (AT d ) 2 
7rr 87„r(l - Td) 
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1 oo — -n 

FIG. 28. Top: Scanning electron micrograph of the which-path device. Bright regions indicate 
metallic gates. A quantum dot is defined in the right arm of a multi-terminal AB ring. A quantum 
point contact to the right of the quantum dot serves as a which-path detector. Bottom: Schematic 
description of the electrodes and contacts in the which-path experiment. The interferometer is 
composed of three different regions, collector C, emitter E, and base regions B. The quantum point 
contact is to the right of the quantum dot. A finite voltage Vd is applied across the quantum point 
contact. Taken from Ref. M. 
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FIG. 29. AB oscillations of the collector current Iq. The solid line is measured with QPC drain 
source voltage Vd = fiV. The dotted line, with reduced visibility, is measured with Vd = 100 fiV. 
Taken from Ref. fli. 
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FIG. 30. Conduction characteristics of the calibration device shown in the inset of (a), (a) 
Conductance through the quantum dot and transmission through the QPC both as a function of 
the plunger gate voltage Vp. (b) Induced average change A7^ in the transmission coefficient of 
the QPC due to adding an electron to the quantum dot as a function of 7^. Each data point is 
obtained by averaging over several Coulomb peaks. Taken from Ref. Q. 
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FIG. 31. (a) QPC transmission coefficient 7^ as a function of QPC gate voltage V g . (b) Mea- 
sured visibility of the AB oscillations as a function of V g for two values of the QPC drain source 
voltage Vd- The visibility is defined as the peak-to-peak signal divided by the average signal. Error 
bars indicate the fluctuations in visibility due to fluctuations of the device's properties. Taken from 
Ref. m. 
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from Ref. [Hi- 
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where T is the intrinsic (zero-temperature) width of the QD resonance. The solid line 
in Fig |3l](b) shows the prediction of Eq. (|5.1|) , where T = 0.5fieV was used as a fitting 
parameter. The quantity ATd was taken from the calibration device, cf. Fig. [J{]. 

The dependence of the visibility on the drain source voltage Vd is illustrated in Fig. |32 . 
For eVd 3> kT, the visibility drops linearly as expected from Eq. (|5.1| ). A deviation from 
the linear dependence occurs near eVd ~ kT. 



2. Experiment of Sprinzak et al. 

The which-path information in the setup of Buks et al. is encoded in the current flowing 
through the QPC-detector. Sprinzak et al. |J used a different approach and devised a 
setup in which the which-path information enters only as a quantum mechanical phase. 
The setup is shown in Fig. [33]: A QPC serves as the detector for the oscillation of charge in 
a double quantum dot (DQD) interferometer. The device is subjected to a high magnetic 
field corresponding to two filled Landau levels of the two-dimensional electron gas. Under 
these conditions most current flows in a chiral motion along the the edges of the sample. 
Scattering between the edge states only occurs at the QPC constriction where the two edges 
of the sample are close together. The nearby DQD interferometer leaves the current through 
the QPC unaffected but does change the phase of the transmitted electron mode. The phase 
of the transmitted electrons therefore contains information about the state of the DQD. Note 
that the DQD is not inserted in a standard two-path interferometer since the presence of 
a high magnetic field would prevent electrons from choosing either path with nearly equal 
probability. The experiment, therefore, does not probe the path of electrons through the 
interferometer but the state of the interferometer. In keeping with Ref. || we still refer to 
the setup which-path experiment. 

The dephasing rate of the interferometer can be extracted from the conductance through 
the DQD. The DQD is being tuned to resonance by means of the two plunger gate voltages 
V p i and V P 2- For the two quantum dots in series, resonances are found when resonant levels 
in both dots are degenerate. The resulting Coulomb peaks are located on a hexagonal lattice 
and each peak has the width 2T where T is the intrinsic resonance width of a single dot 
|6^ ) [12| . In the presence of the detector one expects to observe a reduction of the peak height 



and a broadening of the peak width to 2(T + h/td). A single resonance peak with its contour 
at half maximum is shown in Fig. [34]. Due to the asymmetry of the peak shape, the area in 
this contour is taken as a measure of the dephasing rate. 

Figure shows the dependence of the contour area and the peak height as a function 
of the QPC transmission probability 7^. Both plots were obtained for a fixed QPC voltage 



Vrf.The measured area (Fig. |35|a) yields a dephasing rate proportional T d (l — T d ). Similarly, 



the peak height (Fig. |3lS|b) has an inverse dependence on that expression. Note that the result 
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FIG. 33. Schematic illustration of the DQD interferometer and the QPC detector. The device 
is studied in strong perpendicular magnetic field (5-10 Tesla) that imposes chiral electron motion 
near the edges of the sample. The edge states are partly transmitted (jt)) and partly reflected (|r)) 
by the QPC. The DQD is weakly coupled to its own leads and is tuned to resonance by the two 
plunger gate voltages V p i, V P 2- Taken from Ref. Q. 




FIG. 34. View of one Coulomb peak as a function of the two plunger gate voltages V p \ and 
V P 2- The dashed line is a contour drawn at half maximum of the peak height. The area enclosed 
by this contour is being used as a measure of the peak width. Taken from Ref. || . 
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FIG. 35. Top: The area of the contour at half peak height as a function of the transmis- 
sion probability Td for two values of the applied bias = and 2mV . The dependence agrees 
qualitatively with the expected — Td). Inset: The dependence of the contour area on the 
applied voltage Vd for Td = 0.7. Bottom: The peak height (in units of e 2 /h) as a function of the 
transmission probability Td for two values of applied bias Vd = and 2mV. Taken from Ref. @. 

for the dephasing rate differs from the rate oc (ATd) 2 /[T d [l — T d ] found in the experiment of 
Ref. [||. We show in Sec. [VB| that the difference is due the different nature of the which- 
path information: In the present device this information resides in the transmission phase 
while in the case of the experiment || it is encoded in the transmission probability of the 
QPC. The inset of Fig. ^5| shows that the contour area and hence the dephasing rate grow 
nearly linearly with the detector voltage Vd- 

The measurements show that a phase change in the detector leads to dephasing of the 
interferometer even though no interference experiment is being performed in the detector. 
The question arises onto whether the detector must be phase coherent in order to dephase the 
interferometer. To investigate this question, Sprinzak et al. inserted an artificial dephaser 
(a floating Ohmic contact) in the path of the transmitted electrons before they reached the 
DQD. Even though phase coherence was destroyed, dephasing persisted (however, a small 
decrease of dephasing was observed and traced back to the finite capacitance of the Ohmic 
contact which effectively eliminated high frequency components of the shot noise). 
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B. Theoretical approaches 



In the past, experimental investigations of decoherence (dephasing) have been executed 
with photons, cooled atoms, and neutrons. The which-path experiments of Refs. for 
the first time allowed to investigate controlled dephasing with mesoscopic electron struc- 
tures. Widely different methods have been used [|65| , |66|JT0| ,|8| JTT| j67| H69l^ll to investigate the 
problem theoretically. Conceptually, some of the authors investigated the influence of the 
interferometer on the state of the detector |65|,|8|,|9|] while other authors [^,|l],|Il],|37]] analyzed 
how the detector affects the interferometer. Both descriptions lead to equivalent results as 
is known since the work of Stern et al. 



70]. 



In Sec. |VB 1| we present a simple argument for the dephasing rate in the which-path 
interferometer. The argument displays the intimate link between the dephasing rate and the 
efficiency with which the QPC detector measures the state of the interferometer. Rigorous 
derivations of the dephasing rate and the suppression of AB oscillations in a which-path 
interferometer are given in Sec. [V B 2[-|VB~5 . 



1. Heuristic argument 

A qualitative argument for the dephasing of quantum dot states due to the coupling to 



a QPC detector was given by Aleiner et al. |65[ and Buks et al. ||. Their argument relates 



the dephasing rate to the which-path information obtained in the detector. The original 
argument only considers the information contained in the current through the QPC. Below 
we generalize the argument to include, in addition, the information encoded in the phase of 
the transmission amplitude through the QPC. 

An electron in the quantum dot interacts electrostatically with electrons in the QPC 
detector. The interaction generally modifies both the modulus and the phase of the transmis- 
sion amplitude through the QPC. The change in the transmission probability induced by the 
presence of one electron on the quantum dot is denoted by A 7^; the corresponding change 
of the transmission phase is A(f>. In experiment, these changes cannot be detected with cer- 
tainty due to current and phase fluctuations. The current fluctuations are due to the well- 



known quantum shot noise (for reviews on shot noise see Refs. |7l|-|73|): Let N = (2eVd/h)tq- 
be the number of electrons that probe the QPC biased with the voltage Vd during a time 
t-r- Then the number of transmitted electrons Nq- is a binomial random variable with the 



expectation value (Nq-) = NT d and the standard deviation 5Nq- = J NTd(l — Td). From the 



uncertainty 5Nq- one finds the uncertainty 5T d = y — Xi)/N in 7^. Detection of the 
electron on the quantum dot requires that 5T d is less than or equal to AT d . The equality 
yields the dephasing rate 
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1_ e^JA7^_ 
tr h T d (l - T d ) 

reflecting the which-path information encoded in the detector current. A similar reasoning 
applies to the transmission phase 0: The uncertainty in is given by 5<p ~ where 
SNfj) is the uncertainty in the number of electrons transmitted in the time t^. The condition 
that 5(j) equals the phase change A0 yields the dephasing rate 

Ioc^T d (l-T d )(A0) 2 . (5.3) 

t(f> II 

Adding the dephasing rates due to the change in the current and the phase, one finds that 
the total dephasing rate has the form 

I = c * W +c * S(1 _ Tj)w) , (54) 

t d h TAl - T d ) h 



where Cr and C<f> are coefficients of order 1. The rigorous calculations presented in Sec. V B 
support the simple estimate (|5.4|) and yield expressions for the coefficients Cr and Ca. 



2. Green function 

Controlled dephasing and the suppression of AB oscillations due to the interaction with 



a QPC have been investigated by Aleiner, Wingreen, and Meir | 65fl . The model consists of 
a QD in the Coulomb blockade regime which is embedded in an AB ring and capcitively 
coupled to a QPC (Fig. ^6|). It is assumed that kT is much smaller than the single-particle 



level spacing A. At a transmission peak two charging states of the dot with N and N + 1 
electron are degenerate. Electrons pass through a single level in the dot. The amplitude 
£qd(^) can be expressed in terms of the retarded Green function of this level, 

t QD (e) = -iy/4T L T R J dte' e 'G QD (t), (5.5) 

where Tl,r are the partial width for decay of the level through the left and right tunneling 
barrier, respectively (we put Ti = 1). 

Electrons in the QD interact with electrons in the QPC. This interaction modifies the 
local potential in the QPC. We use the standard description of the QPC as a ID system 
of noninteracting electrons. The Hamiltonian of the QPC when exactly N electrons 
occupy the QD is written in terms of exact scattering states, 

— kWl(k)Mk) + A(k)Mk)}, (5.6) 

where ipL,R are the destruction operators for left and right moving scattering states, respec- 
tively. A summation over spin indices is implied. The presence on the iV + 1 electron changes 
the QPC Hamiltonian to i/jv+i = Hn + V, where 
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QPC 

FIG. 36. Schematic view of the which-path interferometer. The quantum dot is embedded in 
the lower arm of an AB ring. A QPC in the vicinity of the ring acts as a which-path detector. 

V = V LL + V RR + V LR , (5.7) 

with 

Vll(rr) = xj dk ^ 2 ipl {R) {ki)il) L {R) (fa) , (5.8) 

Vlr = XlrJ d -^WU R) {h)^L { R ) {k 2 ) + H.c.]. (5.9) 

The wire with the QPC is connected to two reservoirs between which the finite bias eVd is 
applied. 

The Green function of the resonant level in the dot interacting with the wire is of the 
form 

G QB (t) = -iQ(t)e- iEot - n [P N+1 A4t) + P N A + (t)}, (5.10) 

where E is the energy of the level in the absence of the QPC, and Pn(n+i) the probability 
of the charging state N(N + 1) of the dot, P^ + Pjv+i = 1- The total width is T = Tl + Tr. 
The coherence factors A±(t) describe the response of the wire to the addition (removal) of 
an electron from the dot, 

A + (t) = (5.11) 
A-(t) = (e iHirt e- iH ^*) Hlt+1 . (5.12) 

The expectation values are taken with respect to an equilibrium ensemble in the wire with 
the Hamiltonian and H^ + i, respectively. 

An interpretation of dephasing follows directly from these expressions. They can be 



understood [|70| as the scalar product of two states of the environment (the QPC) at time 



t. These states are identical at t — 0, and thus A±(t = 0) = 1. Then, one state evolves with 
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the Hamiltonian Hn, the other state with i/jv+i- Since the environment is not observed in 
the experiment, its coordinate is being integrated over. Quantum coherence is lost once the 
states of the environment have become almost orthogonal. Then the scalar product becomes 
almost zero and the environment has identified the state of the QD. 

The operators Vll and Vrr mix scattering states that propagate in the same direction. 
This only changes the phase of the transmission amplitude through the QPC The mixing of 
scattering states moving in opposite directions, described by Vlr, modifies the transmission 
coefficient. Note that only the latter contribution is affected by the voltage drop in the wire. 
Hence, the nonequilibrium (finite V d ) part of the dephasing rate is entirely due to this latter 
contribution. 



Aleiner et al. |65[ obtained A± for arbitrary A using known results for the orthogonality 
catastrophe, i.e. the response of an equilibrium noninteracting electron system to a sudden 
perturbation. The nonequilibrium dephasing can be calculated in an expansion to second 
order in \lr- The result has the form 

A + (t) = ( - ^ T Y + \~^Ht, k T,eV d) _ (513) 

\t,o smh irkl 1 1 

Here, £o is an high-energy cutoff, the smaller of the Fermi energy in the wire or the inverse 
rise time of the perturbation in the wire. The definitions of the exponents a, 7 and the 



crossover function h can be found in Ref. [65]. It suffices to note that h(t)/t — > for large 
t. For a symmetric dot, the dephasing rate is given by 

1 = eV* ( A T d f- 

u h mx-r t y ( ' 

which agrees up to a constant factor with the estimate in Eq. ( |5.2|) . 

The calculation of the coherence factor A_(t) yields A_(t) = A + (t)*. At the transmission 
peak, P/v = Pn+i = 1/2, and the transmission amplitude £qd is the Fourier transform ( |5 . 5|) 
of Gqd(£). The final result can be approximated by [pSj] , 



2,r^rir; kT + T„\" kT + r M + eVA' 

( « D = SrT^ [—£-) { 6 j ' (5 ' 15) 

where the total width r to t = T + Td is given in terms of the intrinsic width T and the width 
T d = h/t d introduced by the dephasing. For the which-path experiment || the exponents 
a, 7 appear to be rather small so that the last two terms on the right hand side of Eq. ( |5.15 ) 
yield a factor 1. Moreover, T, <C kT, so that we can expand £qd to find 



Try/rZT^ / eV d (AT d f \ 
" 2kT l, 1 ~ -W32T d (l-T d ) ) ' {5M) 

The term in brackets on the right-hand-side represents the reduction of the visibility due 
to the interaction with the QPC. 
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3. Influence functional 



A conceptually different approach to the which-path interferometer was presented by 
Levinson p6fl . He studied an isolated quantum dot (not connected to leads) with a single 
level at energy Eq and Hamiltonian Hqd = Eqc^c. The QPC and the interaction between 
QPC and QD are modeled precisely as in Sec. |VB2| , so that the total Hamiltonian takes 
the form 

H = H QB + H N + JcV } (5.17) 



where and V are defined in Eqs. ( |5.6|) and ( |5.7|) -( |5T9|) , respectively. Levinson investi- 
gated the effect of the environment on the QD. The charge transmitted through the QPC 
creates a fluctuating potential at the QD that modulates the electron states in the QD. This 
modulation dephases the states in the QD. 

The coherence of the QD is quantitatively described by the expectation value (c(t)) taken 
with respect to the equilibrium ensemble with the Hamiltonian H. Here, c(t) = e tHt ce ~ tHt 
is the amplitude in the Heisenberg representation. The QD has a coherent part if (c(t)) ^ 0, 
while the case (c(t)) = and (ct(t)c(i)) ^ corresponds to a totally incoherent quantum 
dot. The equation of motion for c{t) reads 

*^ = i{E Q + V(t))c{t), (5.18) 
where V(t) = e Ve~ . According to Eq. ( [5.18 ), V{t) may be interpreted as a time 



dependent modulation of the energy Eq. Solving this equation and averaging, one finds 

(c(t)) = (c(0)e- iEot T t e-^o v ^ dt ^, (5.19) 

with the time ordering operator T t . Note that the time ordered exponential on the right 
hand side of this equation is identical with the product exp(iH^t) exp(—iHN + it) that enters 
in the coherence factors A± defined in the previous section. This proves the equivalence 
of the present approach with the one described in Sec. |V B 2|. An expectation value of the 



type of that on the right hand side of Eq. ( 5.19 ) is known as the Feynman- Vernon influence 
functional. 



The average in Eq. ( |5.19| ) is calculated assuming that the process described by 



V(t) is a Gaussian process. This amounts to the approximation (c(t)) = (c(0)) 
x exp(-i£ t)exp[-(l/2)$(t)], where $(t) = f Q t dt'f t dt"K(t' - t") with 

K(t)= 1 -[(V(t)V(0)) + (V(0)V(t)}}. (5.20) 

The correlator K(t) decays in time on a characteristic time scale t^. For large t 3> t^, one 
finds (c(t)) = (c(0)) exp(— iE t) exp(— t/td) where the dephasing rate is given by 
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dtK(t). (5.21) 



The dephasing rate may be found using known results for the current correlation functions 
in a QPC. For a symmetric dot one finds the results 



l/t d ~ AeV d /h for high bias eV d > kT, (5.22) 
l/t d ^ .4(e^) 2 /(/iA;T) for low bias eV d < fcT, (5.23) 



with ^4 = (/S.T d ) 2 /[%ixT d {l — T d )]. The result for the high bias limit is identical with Eq. (|5.14|) 



obtained by Aleiner et al. p5 . 



4- Master equation 

Density matrix approaches to dephasing in coupled quantum dots were devised by 
Gurvitz JTof and Hackenbroich et al. (T^j. The approaches include the dephasing in the 



which-path interferometers as a special case. Gurvitz' work is reviewed in Sec. |VTJ The 
derivation given below is similar in spirit to the approach of Ref . |T| . 

We consider a quantum dot with a single energy level. The QD is coupled to a QPC 
but otherwise isolated from the environment. For the Hilbert space of the QD we choose 
the basis \a), \b), representing an empty and an occupied dot level, respectively. We assume 
that the wire with the QPC supports only a single transverse mode so that the scattering 
matrix S"qpc through the QPC is a 2 x 2 matrix that depends on the occupation of the QD, 

s h a , if the QD is empty, 

I Sb, if the QD is occupied. 

The relations ( |5.24j ) may be combined to the two-particle scattering matrix 

Sera' = &aa> [KaS a + KbSb], (5.25) 

where both a and a' can take the values a, b. Let p tot = p <g> pqpc be the density matrix of 
the total system prior to the passage of an electron through the QPC. We choose pqpc = 
diagfl, 0] representing an incoming particle from one side of the QPC. The density matrix 
after scattering is p' tot = Sp tot S^ . The reduced density matrix p' of the QD is obtained by 
tracing over the QPC variables 

p' = Ti^ c [Sp tot S^. (5.26) 

In the basis of incoming and outgoing scattering states we can parametrize the scattering 
matrix 
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S=\' a " a \ = \ cos ^« ex PN™] ism6 a exp[i(j) t/a \ ^ o _ 

1 ' ' isin6» a exp[i0 ta ] cos # a exp[z0 r / a n 

in terms of one angle 9 a and four phases a . As a consequence of the unitarity of S a , 
the phases are not independent from each other but obey the relation exp[z(0 ra — ta )] = 
— exp[— i(4> r 'a — <fit'a)]- A parameterization similar to Eq. ( p.27| ) but in terms of an angle 



9 b and phases <j) b may be used for the scattering matrix S b . Combining Eqs. ( |5.25| ), (|5.27|) 



and substituting the result into Eq. Q5.26Q one finds that the diagonal elements of p do not 



change upon scattering through the QPC. The off-diagonal elements change according to 

P' ab = V Pah , (5.28) 

where 

V = r* a r b + t* a t b . (5.29) 

Note that p' ba = p'* b since p' is hermitian. 

We will assume that the scattering through the QPC takes place on a time scale much 
shorter than the relevant time scales of the quantum dot. In particular, the scattering time 
shall be much shorter than the decoherence time t d . This assumption allows us to use the 
Markov approximation and neglect memory effects in the QPC. Then the density operator 
after the scattering of n electrons through the QPC is given by 

Pab(t)=V n e^ nAt p ab (0), (5.30) 

where At = h/(2eV d ) is the average time between two scattering events, t = nAt, and 
E at b are the energies of the empty and occupied dot level, respectively. The term involving 
e = (Eb — E a )/h accounts for the free evolution of the QD. Note that V is generally a 
complex number. From its modulus we can read off the dephasing rate 

£ = -imin (5.3D 

while the phase gives rise to the energy shift 

AE = ^argP. (5.32) 

Both the dephasing rate and the energy shift can be obtained explicitly in the weak-coupling 
limit V ~ 1. Then V may be expanded in A# = 9 b — 8 a , A0 r = <f) rb — ra , and A0 t = 
4>tb — <j>ta- Using the parameterization (|5.27|) , we can express A# in terms of the transmission 



coefficients T b;a = \t bA \ 2 . This yields (A#) 2 = (AT d ) 2 / [4T d (l - T d )\ where AT d = T b -T a and 
T d = {% + T a )/2. Substitution into Eqs. (|5.31|) , (|5.32| ) yields the dephasing rate 
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, , rn T , +-^T d (l-T d )(A ( j ) y, (5.33) 
t d h AT d {\ - T d ) h 

and the energy shift 

AE = — *(1 - T d )A0 r + -^T d A0 t , (5.34) 

7T 7T 

where A0 = A0 r — A0 t . The evolution of p in the weak-coupling case may be cast in the 
form of a master equation 



dp 
~dt 



1 AE 

i 



td h 




(5.35) 



The diagonal elements of p do not change in time. Thus the occupation probability of the 
quantum dot does not change in time. By contrast, the off-diagonal elements decay with 
the dephasing rate l/t d . The evolution of p thus indeed reflects dephasing (decoherence) 
and not energy relaxation. 

The dephasing rate ( |5.33| ) reduces to the results derived in Sees. |V B 2[ |V B 3| when the 



phase change A(j) is negligible. The effect of non-zero A0 has been investigated in a paper 
by Stodolsky |B7|| . His result for the dephasing rate (oc 1 — ReV) differs from Eq. ( |5.31| ) and 
does not reproduce the dependence l/t d oc T d {l — T d ) observed experimentally 0. 

5. Charge relaxation 

An new approach to the dephasing problem was developed by Biittiker and coworkers 
75|.[76|j68|,|69| . The authors consider the charge and potential fluctuations that result from the 



Coulomb coupling between mesoscopic conductors. These fluctuations govern the dephasing 
process and thus the dephasing rate. Based on earlier theoretical work on charge relaxation 
in interacting conductors |75|.|T6| , Biittiker and Martin [§8],^9j calculated the dephasing rate 
for a QD coupled to a QPC. The method applies to dephasing both from changes in the 
transmission probability and from changes in the transmission phase. The results obtained 
in both cases coincide with the formulae given earlier. Biittiker and Martin generalized their 
results to the case of more than one channel in the QPC and to a device that includes a 
phase randomizing voltage probe between the QPC and the QD (cf. the experiment of Ref. 
0). It was found that the dephasing rate is unaffected by the voltage probe if there is only 
one channel in the QPC. 
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VI. QUANTUM ZENO EFFECT 



The frequent repetition of a decohering measurement leads to a striking phenomenon 



known as the quantum Zeno effect ]77,7S|: The suppression of transitions between quantum 
states. The standard example is a two-level system with a tunneling transition between the 
two levels. Assume that the system at t = is in prepared in the state a. For small times 
t, the probability to tunnel out of that state is P a _»f,(t) = ^2-|(5|f2o|a)| 2 t 2 , where il /2 is the 
tunneling matrix element, and ujq = Qq/H the tunneling frequency. However, if the interval t 
is split into N sub intervals each with the length t/N, and the system is measured at the end of 
each subinterval, the probability for tunneling is reduced to N P a ^(t / N) = (l/N)P a ^(t). In 
the limit of arbitrary dense measurements, N — ► oo, the system is frozen in the initial state. 
Repeated measurements can thus completely hinder the natural evolution of a quantum 
system. 

Despite considerable theoretical work on the quantum Zeno effect there is only little 
experimental proof for it. An experimental test using an induced hyperfme transition of 
Be ions [79| has been reported. Experiments on optical transitions |80j or atomic Bragg 



scattering have been proposed. Further experimental evidence for the existence of the 
quantum Zeno effect is clearly desirable. Gurvitz [[K|] first pointed out that the quantum 
Zeno effect may be observed in semiconductor microstructures. He studied theoretically the 
tunneling of electric charge between two weakly coupled quantum dots. A QPC located in 
the vicinity of one of the dots served as a non-invasive detector for the charge on that dot 
(see Fig. [371). As expected from the quantum Zeno effect, Gurvitz found that the coupling 
to the QPC suppresses the tunneling oscillations between the two dots. New aspects of the 
problem were identified by Hackenbroich et al. [11]. It was found that an ac voltage in the 
QPC leads to parametric resonance and to a strong reduction of dephasing in the quantum 
dots. The resonance occurs when the frequency ui of the ac signal equals twice the frequency 
ojq of the internal charge oscillations in the double-dot system. The spectral density of the 
detector current was studied by Korotkov p2| , p3| . When the detector is weakly coupled to 
the dots the spectral density has a peak close to the frequency ujq reflecting the oscillation 
of charge between the two quantum dots. For strong coupling the peak disappears and the 
spectral density develops a Lorentzian shape. In this regime charge transfer takes place via 
random jumps rather than by periodic oscillations. 

In Sec. [VI A| we investigate electron tunneling between two quantum dots in the presence 
of a QPC detector. The problem is studied in terms of the reduced density matrix of 



the coupled dots. The spectral density of the detector current is analyzed in Sec. VI B 



In Sec. |V1 C| we modify the system by connecting the dots to external electron reservoirs. 
We show that the dc current through the dots is affected by the detector and provides an 
indirect signature of the quantum Zeno effect. In Sec. [VI D we analyze the decay of an 
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FIG. 37. Mesoscopic device with two quantum dots and a QPC. Electric charge can tunnel 
coherently between the dots. The QPC measures the charge accumulated in one of the dots. 



unstable quantum system, and investigate whether continuous measurements slow down its 
decay rate. 



A. Charge oscillations 

1. Damping 

Consider electric charge oscillating between two quantum dots in the presence of a QPC 
in close vicinity (Fig. |37|). The coupled dots may be viewed as the mesoscopic realization 
of a double-well potential. We assume that there is only one energy level in each dot and 
that the dots are occupied with a single electron. Let the states \a) and \b) with energies E a 
and E b represent the electron in the left and right dot, respectively. We first study the case 
that the interaction with the QPC is negligibly small. Then the dynamics of the two-dot 
system is governed by the tunneling Hamiltonian (Qo/2)(\a)(b\ + \b)(a\). The elements pij 
of the density matrix satisfy the Bloch equations, 

Paa = «y (Pafe ~ Pba) , (6.1) 

Pbb = i-^{Pba - Pab), (6.2) 

Pab=i£pab+i L ^-(Paa- Pbb), (6.3) 

with the frequencies e = (Eb — E a )/h and uq = VLq/U. The diagonal elements p aa (t) and 
Pbb(t), respectively, are the probabilities to find the electron in the left and right dot. The 
off-diagonal elements satisfy p a &(£) = pl a (t) since p is hermitian. Solving the Bloch equations 
for the initial conditions p = diag[l, 0] one finds 

= WHj ^ (6 . 4) 
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where u = yu 2 + e 2 . An electron located initially in the left dot oscillates between the dots 
with the frequency uj. The amplitude of these oscillations is uJq/(uJq + e 2 ). Thus, if e ^> u , 
the electron essentially remains localized in the left dot. 

Now consider the charge oscillations when the QPC detector is active. The transmission 
through the QPC depends on the electron position in the dot. We first assume that the 
electron position only affects the transmission probability and not the transmission phase. 
Let T a (%) be the transmission probability if the electron is in the left (right) dot. Bloch 
equations for this case were derived by Gurvitz |]TD1 . The resulting equations for the diagonal 
elements coincide with Eqs. ( |6.1| ), (|6.2|). Equation ( |6.3| ) is replaced by 

Pab = itPab + ^(Paa ~ Pbb) ~ ^dPab, (6.5) 

with the dephasing rate Kd = eVd(ATd) 2 /[8hTd]- Dephasing causes an exponential damping 
of the nondiagonal matrix element. Combination of Eqs. (|6.1|) , ( |6.2| ), and ( |6.5| ) determines 
the time evolution of p. For weak damping ujq > Hd and e = 0, one finds damped oscillations 



Paa(t) - 1/2 ~ exp[-K d t] cos J u 2 - K 2 d t. (6.6) 



Note that the damping causes a redshift of the frequency J oj 2 — k 2 , and, hence, slows down 
the charge oscillations between the two dots. This suppression of transitions due to measur- 
ing with the QPC is a manifestation of the quantum Zeno effect. The off-diagonal elements 
of p vanish for large times in agreement with the standard picture of decoherence, and the 
density matrix reduces to the random statistical ensemble p — > diag[l/2, 1/2]. 

Gurvitz' derivation of the rate equation implicitly assumes <C 1. The generalization 



to arbitrary value < Td < 1 was obtained by Hackenbroich et al. \\ 1] using a two-particle 



scattering approach similar to the calculation of Sec. |V B 4|. The resulting Master equation 



is identical with Eqs. (|6.1|), (|6.2| ), (|6~5| ) but yields the dephasing rate 



eV d (AT,) 2 

K d = 6.7 

h 8T d (l - T d ) 

that differs from Gurvitz' result by a factor 1/(1 — T d ). The result (|6.7|) holds provided the 
quantum dots only affect the transmission probability through the QPC. The result may 
easily be generalized to include a change in the transmission phase through the detector. 
One then finds K d = l/(2t d ) where l/t d is the dephasing rate ( |5.33D of a single quantum dot 
computed in Sec. |VB| . [We note that the problem considered here reduces to the dephasing 
in a single quantum dot in the limiting case u>o = 0.] 

2. Parametric resonance 
The charge oscillations between the two quantum dots may display a parametric reso- 



nance if the voltage drop Vd has an ac-component Consider Vd(t) — Vq — V\ sin toit with 
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FIG. 38. The occupation probability of one quantum dot as a function of time. The curves 
correspond to different values of the dephasing rate: (dashed), = 4loo (dot-dashed), and 
Kd = lGu>o (solid). The figure is obtained for e = 0. Taken from Ref. [p7cf| . 




FIG. 39. Parametric resonance in the double-dot system coupled to a QPC-detector. The upper 
part shows the oscillations of p aa for constant voltage (dashed curve) and for a time-dependent 
voltage as shown in the lower part (full curve). Damping of the oscillations is reduced by a factor 



two. Taken from Ref. [11|. 



Vq, V\ > and V\ < Vq. Substitution into the equation of motion for p aa yields a damped 
harmonic oscillation with an oscillatory damping constant. Using a simple ansatz for p aa {t) 
and neglecting terms of order 0{V^), one is led to an equation of the Mathieu type which is 
known [5S| to display parametric resonance close to the frequencies ui\ = 2u/n where n is a 
positive integer. Parametric resonance is most pronounced for uj\ ^ 2u. The damping near 
the resonance is strongly reduced, Kd = (A7^) 2 /[8h%i(l — %)]e{Vo — Vi/2). The resulting 
time evolution of p aa near resonance is illustrated in Fig. [39] and compared with the case 
where is time-independent. 



B. Detector current 

While the reduced density matrix p of the two quantum dots is a convenient tool for 
the theoretical investigation of the quantum Zeno effect, one cannot directly measure p 
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in experiments. However, information about the dynamics in the quantum dots may be 
obtained from the current passing through the detector. The power spectrum of the detector 
current was studied by Hackenbroich et al. [O] and by Korotkov \ 82 ,l 



Korotkov's approach p2| , p3| is based on the description of decoherence in terms of a 
stochastic wave function. This approach is well-known in quantum optics ||84|| . Variants of 
the approach are known as quantum trajectory or quantum jump approach. The method 
describes the evolution of the quantum system under investigation (the two quantum dots) 
conditioned on a particular measurement of the detector system (the QPC). For the condi- 
tional evolution of the quantum dot density matrix given the detector current I(t) one finds 



82.83 



Paa = ^(Pab ~ Pba) ~ ^r" PaaPbb[l '(*) ~ h] , (6.8) 
2 D 

Pab = itPab + ^{Paa ~ Pbb) + -^(Paa ~ Pbb)[I(t) - I Q ]p ab - ^p ab , (6.9) 
I D 

where S is the low-frequency spectral power of the detector shot noise, and J = (^a+A)/2, 
AI = lb — I a where I a and h are the average currents when the electron localized on the left 
and right dot, respectively. The damping rate 7 accounts for dephasing due to interaction 
with some environment not included in the detector. Eqs. ( |6.8|) , (|6.9|) are supplemented by 
the equation 

I(t) -I = AI(p aa - p bb )/2 + f (t) (6.10) 

for the detector current. Here, £(i) describes a Gaussian random process with zero average 
and spectral density 6V Equations (|6.8|)-( |67L0"|) describe the evolution of the system for a 



single realization of the detector current. Averaging over £ reduces these equations to the 



deterministic Master equation for the evolution of the ensemble average derived in Sec. VIA . 

The evolution of I(t) and p aa are illustrated in Fig. ^D] for two different values of the 
coupling a = Kd/^o- Note that the results hold for a particular realization of the random 
process £(t). Figure (5]( a ) is obtained for a = 0.1 representing a weakly coupled detector. 
Charge can oscillate several times between the two dots while generating only weak modu- 
lations of the detector current. Superimposed on the modulations is the detector shot noise. 
Fig. |40|(a) shows the detector current averaged over different time intervals. For small in- 
tervals the signal is noisy, while for long intervals individual oscillations cannot be resolved. 
As a result, charge oscillations are difficult to observe in the weak-coupling limit. The situ- 
ation is different in the strong coupling case a > 1 (Fig. 0(b)). The strong influence of the 
detector suppresses the quantum oscillations, so that the charge performs random jumps 
between two localized states. In this case, the properly averaged detector current follows 
the evolution of the quantum dots. 
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FIG. 40. Evolution of p aa (thick line) and the corresponding detector current I(t) (thin slid, dot- 
ted, dashed line) averaged using rectangular time windows of size r a . (a) Weak-coupling a = 0.1, 
(b) strong-coupling a = 5.0. Taken from Ref. 
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The detector current may be further characterized by its spectral power. For a symmetric 
double-dot (e = 0), the spectral power can be obtained analytically |83| yielding 



2 4 

S{u) = S + 4S-0— 2 2 2 . (6.11) 

Figure f|I] shows S^u;) for different values of the coupling strength a. For weak coupling, 
the spectral power displays a peak close to the frequency ujq. The peak height equals four 
times the noise background. The full width at half maximum is given by 2n d and the peak 



is centered at the shifted frequency ytUg — K d [|TT|]. The frequency shift with respect to the 
frequency ujq is a consequence of the quantum Zeno effect. The peak in S(u) gradually 
disappears with increasing coupling. In the strong coupling limit, S(u) has a Lorentzian 
shape known from the classical theory of telegraph noise |53] . 



C. Current through the quantum dots 

Due to high-frequency noise, a frequency resolved measurement of the QPC current is 
difficult to perform. To avoid the problems of high-frequency measurements several authors 
T0| , [TT|| investigated manifestations of the quantum Zeno effect in dc transport. The dc 



current through two quantum dots with a QPC in the vicinity was studied both for the 
dots in series and with the dots in a parallel circuit. The current displays a characteristic 
dependence on the detector efficiency which provides an indirect signature of the quantum 
Zeno effect. 



1. Dots in a series 

Consider the transmission through a double-well structure consisting of two coupled 
quantum dot in a series. The dots are modeled as in Sec. |V1 A| in terms of two levels with 
energies E a , E b , and coupling matrix element Qq/2. Dot a is coupled to an electron reservoir 
a and dot b to a reservoir b. The coupling strength is described by the partial widths T a ^ 
for the decay of level a(b) into the reservoir a(b). The Fermi energies in the reservoirs are 
chosen such that E b F <C E a ,E b <C E F . The current I s through the double-dot system can 
be calculated fn| using rate equations and is given by 



(r 6 + TiK d )VL 



ft A{E a - E b y + (r„ + hn d y + n 2 (r b + n Kd ) (a + ± 



(6.12) 



For small energy difference \E a — E b \ the current decreases with n d . However, for large 
\E a — E b \ the current increases with n d . This may be understood QTUJ as resulting from the 
derealization of the electron due to a continuous measurement of the charge in one of the 
dots.. 
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FIG. 42. Two coupled quantum dots form a double well potential and interact with a QPC 
measuring the position of an excess electron in the two dots. A current is driven through the dots 
via four external leads. The transmission and reflection coefficients of the leads are denoted by T u , 
T h Ri. 

2. Parallel quantum dots 



We now discuss the quantum Zeno effect in a set-up for dc transport across two parallel 
quantum dots. The arrangement, proposed by Hackenbroich et al. ]TT[, is shown in Fig. |4"2| . 
Two dots are coupled via a tunneling barrier. Each dot is connected with two external leads 
that allow for the transport of current to and out of the dot. With current flowing into the 
lower dot, we calculate the branching ratio of the current transmitted through the upper 
and the lower dot, respectively. 

We assume that both dots are in the Coulomb blockade regime, and that kT is much 
smaller than the single-particle level spacing in the dots. Then we need to consider only a 
single energy level in each dot. Both levels are degenerate with energy E — i/2T. The width 
T is due to the coupling to the leads. A QPC measures the charge in the lower quantum dot. 
The measurement modifies the transmission through the two-dot system. This modification 
is found by calculating the two-particle scattering matrix for an electron passing through 
the QPC and another electron oscillating in the two-dot system. 

We model the QPC in terms of plane waves with energy e^, mean density pp, and 
Hamiltonian 



#qpc + V= 5> fc &t& fc + Vwdjdjblbv . 



(6.13) 



k,k' 



Here, V describes the interaction between the QPC and the lower QD. The potential Wkk 1 
mimics the constriction in the wire with the QPC while d\ (b\) denotes the creation operator 
for the state on the lower QD (for the QPC plane wave states). 

In each of the leads, we consider only one transverse channel labeled by c = a, p, with 
a = l,u denoting the lower (1) and upper (u) lead, and // = +,— the left and right lead. 
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The two-particle scattering matrix can be obtained from the Lippmann-Schwinger equation 
yielding the scattering amplitudes 



S C c',kk' = Scc'bkk' - ^ijclc'Gaa'M'- (6.14) 

Here, 7 C is the matrix element for tunneling from channel c to the adjacent dot. The two- 
particle Green function G for the joint transition between dots and QPC is given in terms 
of its inverse 

f/ ~i-u [G^dkk' — Wkk' ^o/2 5kk' \ 

(L* )au',kk'-\ , x • (b.lOj 

y i'o/ 2 °kk i Vkk' — vy kk > — v k k'j 

On the right hand side, we have explicitly displayed the matrix in cr-space. Here Gq 1 = 
E — Eq — e k + iT/2 is the inverse propagator of the single Breit-Wigner resonance, and 
E is the sum energy of the two incoming particles. We have assumed that all 7 C = 7 are 
identical and used the relation T = 4^77*. For V = 0, we can diagonalize the matrix 
tk^kk' + Wkk' by a unitary transformation in fc-space, and G reduces to the product of the 
unit matrix in fc-space and the two fc-dependent coupled Breit-Wigner resonances for the 
double-dot system. Then, all scattering processes are elastic, and the branching ratio for 
the transmission through the upper and lower lead is T^/T^ = tt 2 /T 2 . 

The full two-particle scattering matrix ( |6.14|) allows for energy exchange between dots 



and detector: In contrast to the sum energy of the two incoming particles, the energy of 
electrons in the QPC is not conserved in the scattering process. Such inelastic processes are 
essential to ensure the unitarity of the S-matrix. Physically, the energy exchange allows for 
a position measurement of the dot electron without violation of the Heisenberg uncertainty 
relation. 

To calculate the transmission and reflection coefficients through the double-dot system, 
we restrict ourselves to constant scattering potentials Wkk' = W and Vkk' = V . We expand 
G in Eq. ( |6.15| ) in powers of V and resum the resulting series. We obtain two contributions 
to G. The first is independent of V and describes independent elastic scattering through the 
QPC and the dots. This term obviously does not contribute to the quantum Zeno effect. 
The second contribution G involves energy exchange Q = e k — between dots and QPC. 
For fixed incident energy E = E + e k , 

<W=^(«) I . no(2 n + ir) ir (2Si + ir)) ' < 616 > 



with the amplitude 



A ^ Q) = ~f w f w+v ^ 2 + n§)((2fi + ivy - n§) • (6 - 17) 
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Here, Fw = 1 + 2iriWpF. The transmission and reflection coefficients are calculated by 
adding the two-particle scattering probabilities and tracing over the degrees of freedom of the 
QPC. The resulting expressions are further simplified by using a weak-coupling expansion 
to second order in V. In this limit the application of a drain source voltage across the 
QPC is equivalent to the simultaneous scattering of ItVdPF particles in different longitudinal 
QPC modes. The total effect of these particles is obtained by multiplying the result for one 
QPC-particle with the number of longitudinal modes. 

The ^-dependent correction to the branching ratio arises both from coherent (elastic) 
and incoherent (inelastic) scattering. We find that measurements with the QPC-detector 
have a twofold effect: (i) They suppress tunneling from the feeding lead into the lower 
dot and (ii) they suppress tunneling from the lower into the upper dot. Observation (i) is 
reflected in the increase of the reflection coefficient, while (ii) follows from the decrease of 
the branching ratio 



T u r,(o) 



% r/ 0) 



eV d (AT,) 2 
ttT 4T d (l - T d ) 



(6.18) 



-i 

Here, the ratio T^/Tj was given above, and AT, and % were defined in Sec. |V|. Both 
effects (i), (ii) have an obvious interpretation as manifestations of the quantum Zeno effect. 
The second term in the square bracket is up to a factor Y/{2%) equal to the dephasing rate 
found for the isolated double-dot system. 



D. Decay rate 



In early work on the quantum Zeno effect it was argued [FT]] that an unstable quantum 
system slows down its decay under frequent or continuous observations. However, despite 
of further theoretical work |86|] the suppression of decay remained controversial. Indeed, 



a simple argument suggests that the decay rate should not be influenced by observation: 
Consider the exponential decay from a discrete initial state into a continuum of final states. 
The probability for decay during a small time t is linear in t and given by P(t) = Tt/h. For 
N measurements the decay probability is NP(t/N) = P(t) and hence identical to the decay 
probability in the absence of any measurement. 

Elattari and Gurvitz [|8"T| , |8lf| addressed the problem for the system shown in Fig. A 
quantum dot with a single energy level Eq is coupled to an electron reservoir. The states E a 
in the reservoir are dense. The dot is placed near a QPC connected to two separate reservoirs 
filled up to the Fermi energies p^ and p^, respectively. The transmission probability through 
the QPC depends on the occupation of the dot level. Thus, the QPC continuously monitors 
the dot occupation. The dynamics of the entire system is determined by the Schrodinger 
equation ihdt\^f(t)) = H\^(t)), where the Hamiltonian H comprises the quantum dot, the 



S3 



QPC, their mutual interaction and the electron reservoirs. As the initial state |^(0)) we 
choose an occupied quantum dot level and electron reservoirs filled up the Fermi levels /il 
and jiR. 

The time evolution of the system can be expressed in the Bloch-type equations p? 



r 

o"oo = "^oo, (6.19) 

ff«o = i^-^Oa - CTao), (6-20) 

Eq — E a tt a T + hn d 

cr a o = i z CF a Q-i—a 00 — — a a0 . (b.Zl) 

n n 2n 

Here, <t o is the probability for the electron to occupy the dot level, a aa is the probability 
for a transition into level a in the continuum, and Kd is the decoherence rate induced by the 
detector. The coupling matrix element between the dot level and the level a is denoted by 
Q a , and V = 2ir\Q a \ 2 pis the total width for decay into the continuum, where p is the density 
of states in the reservoir coupled to the quantum dot. From Eq. ( |6.19| ) one immediately finds 

a 00 (t) = exp(-rt/H). (6.22) 

Hence, continuous monitoring of the unstable system does not slow down its exponential 
decay. In contrast, the energy distribution of the electron in the continuum, P(E a ) = 
&aa(t — > oo) is affected by the detector. From the solution of Eqs. (|6.19| )- (|6T2T| ) in the limit 
t — > oo, one finds a Lorentzian distribution centered around E a = E , 

The measurement results in a broadening of the linewidth from V to T + hna- The broadening 
does not affect the decay rate V. 
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FIG. 43. Schematic illustration of a QPC-detector near a quantum dot that is coupled to a 
reservoir. The energy level Eq of the dot is occupied by an electron which can decay into the states 
E n of the reservoir. Taken from Ref. 



VII. CONCLUSION 

Semiconductor quantum dots allow to measure and control current at the single-electron 
level. Multiple quantum dots, wires and quantum point contacts can be integrated on a 
single chip into more complicated structures. Such integrated microstructures not only offer 
control over the number of transmitted electrons, they also allow to measure and manipulate 
the quantum state of these electrons. 

We demonstrated that devices employing a quantum dot in an Aharonov-Bohm ring 
allow to measure the phase of the transmission amplitude through the quantum dot. Con- 
trolled dephasing of quantum dot states can be achieved by means of a quantum point 
contact in close proximity to the dot. It is a formidable task to extend these experiments 
to quantum dots in the Kondo regime. Theoretical predictions for the phase shift in this 
regime have been made ||89|| . Novel experiments testing the dephasing of charge tunneling 



between coupled quantum dots have been proposed [[K],|ri[]. The research in this field al- 
lows for the first time to test fundamental principles of quantum mechanics, such as the 
particle-wave duality, with solid state devices. Questions that previously could only be 
addressed in atomic physics or quantum optics may now be amenable |9(J to solid state 
physics. Further question may be investigated, e.g. exploring the connection of dephasing 



and the Fermi statistics [PTJ. These considerations promise interesting future research on 



electron coherence in a wide variety of mesoscopic systems. 
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